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Abstract. We consider the natural Langevin dynamics which is reversible with respect 
to the mean-field plane rotator (or classical spin XY) measure. It is well known that this 
model exhibits a phase transition at a critical value of the interaction strength parame- 
ter K, in the limit of the number N of rotators going to infinity. A Fokker-Planck PDE 
captures the evolution of the empirical measure of the system as N — > oo, at least for 
finite times and when the empirical measure of the system at time zero satisfies a law of 
large numbers. The phase transition is reflected in the fact that the PDE for K above 
the critical value has several stationary solutions, notably a stable manifold - in fact, 
a circle - of stationary solutions that are equivalent up to rotations. These stationary 
solutions are actually unimodal densities parametrized by the position of their maximum 
(the synchronization phase or center). We characterize the dynamics on times of order 
N and we show substantial deviations from the behavior of the solutions of the PDE. In 
fact, if the empirical measure at time zero converges as N — > oo to a probability measure 
(which is away from a thin set that we characterize) and if time is speeded up by iV, 
the empirical measure reaches almost instantaneously a small neighborhood of the stable 
manifold, to which it then sticks and on which a non-trivial random dynamics takes 
place. In fact the synchronization center performs a Brownian motion with a diffusion 
coefficient that we compute. Our approach therefore provides, for one of the basic sta- 
tistical mechanics systems with continuum symmetry, a detailed characterization of the 
macroscopic deviations from the large scale limit - or law of large numbers - due to finite 
size effects. But the interest for this model goes beyond statistical mechanics, since it 
plays a central role in a variety of scientific domains in which one aims at understanding 
synchronization phenomena. 
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1. Introduction 

1.1. Overview. In a variety of instances partial differential equations are a faithful ap- 
proximation - in fact, a law of large numbers - for particle systems in suitable limits. This 
is notably the case for stochastic interacting particle systems, for which the mathematical 
theory has gone very far [24] . The closeness between the particle system and PDE is 
typically proven in the limit of systems with a large number N of particles or for infinite 
systems under a space rescaling involving a large parameter N - for example a spin or 
particle system on % d and the lattice spacing scaled down to jt - and up to a time horizon 
which may depend on N. Of course the question of capturing the finite N corrections has 
been taken up too, and the related central limit theorems as well as large deviations prin- 
ciples have been established (see [23] and references therein). Sizable deviations from the 
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law of large numbers, not just small fluctuations or rare events, can be observed beyond 
the time horizon for which the PDE behavior has been established and these phenomena 
can be very relevant. 

The first examples that come to mind are the ones in which the PDE has multiple 
isolated stable stationary points: metastability phenomena happens on exponentially long 
time scales [29]. Deviations on substantially shorter time scales can also take place and 
this is the case for example of the noise induced escape from stationary unstable solutions, 
which is particularly relevant in plenty of situations: for example for the model in [301 
Ch. 5] phase segregation originates from homogeneous initial data via this mechanism, 
on times proportional to the logarithm of the size of the system. The logarithmic factor 
is directly tied to the exponential instability of the stationary solution (see [30] for more 
literature on this phenomenon). Of course, the type of phenomena happen also in finite 
dimensional random dynamical systems, in the limit of small noise, but we restrict this 
quick discussion to infinite dimensional models and PDEs. 

In the case on which we focus the deviations also happen on time scales substantially 
shorter than the exponential ones, but the mechanism of the phenomenon does not involve 
exponential instabilities. In the system we consider there are multiple stationary solutions, 
but they are not (or, at least, not all) isolated, and hence they are not stable in the standard 
sense. Deviations from the PDE behavior happen as a direct result of the cumulative effect 
of the fluctuations. More precisely, this phenomenon is due to the presence of whole stable 
manifold of stationary solutions: the deterministic limit dynamics has no dumping effect 
along the tangential direction to the manifold so, for the finite size system, the weak noise 
does have a macroscopic effect on a suitable time scale that depends on how large the 
system is. We review the mathematical literature on this type of phenomena in § 11.61 
after stating our results. 

Apart for the general interest on deviations from the PDE behavior, the model we 
consider - mean-field plane rotators - is a fundamental one in mathematical physics and, 
more generally, it is the basic model for synchronization phenomena. Our results provide 
a sharp description of the long time dynamics of this model for general initial data. 

1.2. The model. Consider the set of ordinary stochastic differential equations 



with j = 1,2, . . . , N , {Wj"}j=i,2,... is an HD collection of standard Brownian motions and 



give origin to a diffusion process on S N , where § is the circle R/(2ttZ). 

The choice of the interaction potential </(•) is such that the (unique) invariant proba- 
bility of the system is 





i=l 






where A at is the uniform probability measure on S . Moreover, the evolution is reversible 
with respect to itn Ki which is the well known Gibbs measure associated to mean- field 
plane rotators (or classical XY model). 
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We are therefore considering the simplest Langevin dynamics of mean-field plane rota- 
tors and it is well known that such a model exhibits a phase transition, for K > K c := 1, 
that breaks the continuum symmetry of the model (for a detailed mathematical physics 
literature we refer to [6]). The continuum symmetry of the model is evident both in the 
dynamics (jl.ip and in the equilibrium measure fjl .2[) : if }t>o,j=i,...,N solves (jl.ip . so 
does {(pi' + c}t>oj=i,...,jVi c an arbitrary constant, and ttn,k®c 1 = ^n^k, where G c is 
the rotation by an angle c, that is (Q c (p)j = ipj + c for every j. 

1.3. The N — > oo dynamics and the stationary states. The phase transition can be 
understood also taking a dynamical standpoint. Given the mean-field set up it turns out 
to be particularly convenient to consider the empirical measure 

1 - 

W,t(M) := ^E^W' (1-3) 
i=i 

which is a probability on (the Borel subsets of) S. It is well known, see [6] (for detailed 
treatment and original references), that if /xjvo converges weakly for ./V — > oo, then so 
does fiN,t for every t > 0. Actually, the process itself t \-t {/ijv,t}, seen as an element of 
C°([0, T], where T > and Ai± is the space of probability measures on § equipped 

with the weak topology, converges to a non-random limit which is the process that con- 
centrates on the unique solution of the non-local PDE (* denotes the convolution) 

dtPt(0) = \d 2 e p t (8)-d e {(J* Pt )(0)p t (0)), (1.4) 

with initial condition prescribed by the limit of {hn,o}n=i,2,...- If such a limit probability 
does not have a (C 2 ) density (with respect to the uniform measure), one has to interpret 
(|1.4|) in a weak sense, but actually, even if the initial datum is just in .Mi, that is if it 
does not admit a density or if such a density is not smooth, the probability measure that 
solves (ll.4j) has a density pt(-) £ C°° for every t > 0, see [21] . We insist on the fact that 
Pt(-) is a probability density: f s pt(0) d8 = 1. We will often commit the abuse of notation 
of writing p{9) when p £ Mi and p has a density. Much in the same way, if p(-) is a 
probability density, p, or p(d6), is the probability measure. 

It is worthwhile to point out that (J * p)(0) = — ?R(pi)K sin(#) + Q(pi)K cos(#) with 
pi := f s p(0) exp(i6) d6. This is to say that the nonlinearity enters only through the first 
Fourier coefficient of the solution, a peculiarity that allows to go rather far in the analysis 
of the model. Notably, starting from this observation one can easily (once again details 
and references are given in [6]) see that all the stationary solutions to (II. 4h . in the class 
of probability densities, can be written, up to a rotation, as 

exp(2Kr cos(6Q 
q{9) - = 2nI Q (2Kr) ' (L6 > 

where 2TrlQ(2Kr) is the normalization constant written in terms of the modified Bcssel 
function of order zero (Ij(x) = (2ir)~ 1 J§(cos 0) 3 1 exp(x cos(#)) d6>, for j = 0,1) and r is a 
non-negative solution of the fixed point equation r = fy(2Kr), with ^(x) = I\(x)/ Iq(x). 
Since *&(•) : [0, oo) — > [0, 1) is increasing, concave, ^(0) = and ^'(0) = 1/2 we readily see 
that if (and only if) K > 1 there exists a non-trivial (i.e. non-constant) solution to (|1.4p . 
Let us not forget however that ^(0) =0 implies that r = is a solution and therefore the 
constant density 5- is a solution no matter what the value of K is. From now on we set 
K > 1 and choose r = r(K), the unique positive solution of the fixed point equation, so 
that the probability density q(-) in (|1.5j) is non trivial and it achieves the unique maximum 
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Figure 1. The evolution limit evolution (|1.4[1 instantaneously smoothens an arbitrary 
initial probability and, unless the Fourier decomposition such an initial condition has 
zero coefficients corresponding to the first harmonics (the hyperplane U), it drives it to 
a point Poo - a synchronized profile - on the invariant manifold M and of course it stays 
there for all times. This has been proven in |21] , here we are interested in what happens 
for the finite size - N - system and we show that the PDE approximation is faithful up 
to times much shorter than N: on times proportional to N synchronization is kept and 
the center of synchronization ip performs a Brownian motion on S. 

at and the minimum at tt. Note that the rotation invariance of the system immediately 
yields that there is a whole family of stationary solution: 

M = {<fy(.) : := q{- - ifj) and ^e§}, (1.6) 

and, when qx(-) of course means q X mod(2ir)(')- M, which is more practically viewed 

as a manifold (in a suitable function space, see § 12.21 below), is invariant and stable for the 
evolution. The proper notion of stability is given in the context of normally hyperbolic 
manifolds (see [32] and references therein), but the full power of such a concept is not 
needed for the remainder. Nevertheless let us stress that in [21] one can find a complete 
analysis of the global dynamic phase diagram, notably the fact that unless po(-) belongs 
to the stable manifold U of the unstable solution J- - the solution corresponding to r = 
in (jl.5p - pt(-) converges (also in strong norms, controlling all the derivatives) to one of 
the points in M, see Figure [TJ There is actually an explicit characterization of U: 

U = jp G Mi : J exp(i9)p(d9) = oj . (1.7) 

As a matter of fact, it is easy to realize that if po(-) £ U then (|1.4|) reduces to the heat 
equation dtpt(0) = \d^pt(9) which of course relaxes to 

1.4. Random dynamics on M: the main result. In spite of the stability of M, q^(-) 
itself is not stable, simply because if we start nearby, say from q^>, the solution of fjl .4|) 
does not converge to gy>(-)- The important point here is that the linearized evolution 
operator around q(-) £ M (q is an arbitrary element of M, not necessarily the one in (|1.5|) : 
the phase ip of q^ is explicit only when its absence may be misleading) 

L q u{9) := -u" -[uJ*q + qJ* u}' , (1.8) 

with domain {u 6 C 2 (S,M) : J g n = 0} is symmetric in H_i i/ q - a weighted H-\ Hilbert 
space that we introduce in detail in Section [2.11 - and it has compact resolvent. Moreover 
the spectrum of L q , which is of course discrete, lies in (— oo,0] and the eigenvalue has 
a one dimensional eigenspace, generated by q' . So q' is the only neutral direction and it 
corresponds precisely to the tangent space of M at q(-): all other directions, in function 
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space, are contracted by the linear evolution and the nonlinear part of the evolution does 
not alter substantially this fact |214 [25] . 

Let us now step back and recall that our main concern is with the behavior of 
with N large but finite, and not (ll.4p . In a sense the finite size, i.e. finite N, system is 
close to a suitable stochastic perturbation of f 1 1 . 4 [> : the type of stochastic PDE, with noise 
vanishing as N — > oo, needs to be carefully guessed pj5], keeping in particular in mind 
that we are dealing with a system with one conservation law. We will tackle directly (II. ip . 
but the heuristic picture that one obtains by thinking of an SPDE with vanishing noise 
is of help. In fact the considerations we have just made on L q suggest that if one starts 
the SPDE on M, the solution keeps very close to M, since the deterministic part of the 
dynamics is contractive in the orthogonal directions to M, but a (slow, since the noise is 
small) random motion on M arises because in the tangential direction the deterministic 
part of the dynamics is neutral. This is indeed what happens for the model we consider for 
N large. The difficulty that arises in dealing with the interacting diffusion system (11. lj) is 
that one has to work with (|1.3|) . which is not a function. Of course one can mollify it, but 
the evolution is naturally written and, to a certain extent, closed in terms of the empirical 
measure, and we do not believe that any significative simplification arises in proving our 
main statement for a mollified version. Working with the empirical measure imposes a 
clarification from now: as we explain in Section [2.11 and Appendix [A] if [i and v 6 Mi, 
then fj, — v can be seen as an element of H-i (or, as a matter of fact, also as an element 
of a weighted space). 

Here is the main result that we prove (recall that K > 1): 

Theorem 1.1. Choose a positive constant tj and a probability po G Mi \ U . If for every 
e > 



lim P | \ m -po _i<£ = 1, (1-9) 

then there exist a constant ip$ that depends only on po(-) and, for every N , a continuous 
process {Wn, t }t>o, adapted to the natural filtration of {W ] N }j = i t 2,...,N , such that Wjv,- £ 
C°([0, Tf]; M) converges weakly to a standard Brownian motion and for every e > 



lim P SUp \\fl N ,rN ~ qip +D K w NT \\i < £ = 1) (1-10) 

N -*°° V re [ £ ~< r /] ' / 

where en '■= C/N , C = C(K,po,e) > 0, and 



D K := —j= 1 (1.11) 

'l - (WKr))- 2 

The result is saying that, unless one starts on the stable manifold of the unstable solution 
(see Remark [23] for what one expects if j>o £ U), the empirical measure reaches very quickly 
a small neighborhood of the manifold M: this happens on a time scale of order one, as 
a consequence of the properties of the deterministic evolution law (jl.4p (Figured]), and, 
since we are looking at times of order N, this happens almost instantaneously. Actually, 
in spite of the fact that the result just addresses the limit of the empirical measure, the 
drift along M is due to fluctuations: the noise pushes the empirical measure away from 
M but the deterministic part of the dynamics projects back the trajectory to M and the 
net effect of the noise is a random shift - in fact, a rotation - along the manifold (this is 
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taken up in more detail in the next section, where we give a complete heuristic version of 
the proof of Theorem II. IB . 

Remark 1.2. Without much effort, one can upgrade this result to much longer times: 
if we set Tf(N) = N a with an arbitrary a > 1, there exists an adapted process T 
converging to a standard Brownian motion such that 



lim P sup 



7V-s>oo 



VN,tN°- — Qip +D K N a - 1 W^ [ T 



<£ I =1. (1.12) 



This is due to the fact that our estimates ultimately rely on moment estimates, cf. Sec- 
tion [3 These estimates are obtained for arbitrary moments and we choose the moment 
sufficiently large to get uniformity for times O(N), but working for times 0(N a ) would 
just require choosing larger moments. We have preferred to focus on the case a = 1 this is 
the natural scale, that is the scale in which the center of the probability density converges 
to a Brownian motion and not to an "accelerated" Brownian motion (this is really due to 
the fact that we work on E> and marks a difference with [9j H] where one can rescale the 
space variable). 

1.5. The synchronization phenomena viewpoint. The model (jl.ip we consider is 
actually a particular case of the Kuramoto synchronization model (the full Kuramoto 
model includes quenched disorder in terms of random constant speeds for the rotators, 
see [H [6] and references therein). The mathematical physics literature and the more 
bio-physically oriented literature use somewhat different notations reflecting a slightly 
different viewpoint. In the synchronization literature one introduces the synchronization 
degree r^,t and the synchronization center *jv,t via 

r Njt exp(i& Ntt ) := -^^exp(^ ,7V ) ^= exp(i0)//jv,t( d0)j , (1.13) 

which clearly correspond to the parameters r and ip that appear in the definition of M, 
but rjv,t and *&N,t are defined for N finite and also far from M. Note that if fjl ,9|) holds, 
then both rjv,t and *&N,t converge in probability as N — > oo to the limits r and if), with 
rexp(i?/>) = L exp(i#)po( d#) and the assumption that po U just means r ^ 0. Here is a 
straightforward consequence of Theorem 11.11 

Corollary 1.3. Under the same hypotheses and definitions as in Theorem \1.1\ we have 
that the stochastic process *&n,n- £ C°([e, rj]; S) converges weakly, for every e G (0, Tf], to 
(V>o + D K W.)mod(2ir). 

It is tempting to prove such a result by looking directly at the evolution of ^N,t- 
d*7V,t = ( -K + 



1 \ 1 N 



1 N 

+ -J>s(^-*^)d^)- (1.14) 

But this clearly requires a control of the evolution of the empirical measure, so it does 
not seem that (|1.14|) could provide an alternative way to many of the estimates that we 
develop, namely convergence to a neighborhood of M and persistence of the proximity 
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to M (see Section and Section On the other hand, it seems plausible that one 
could use (|1.14fl to develop an alternative approach to the dynamics on M, that is an 
alternative to Section HI While this can be interesting in its own right, since the notion 
of synchronization center that we use in the proof and *f?Nt are almost identical (where 
they are both defined, that is close to M) we do not expect substantial simplifications. 

1.6. A look at the literature and perspectives. Results related to our work have been 
obtained in the context of SPDE models with vanishing noise. In j8j [15] one dimensional 
stochastic reaction diffusion equations with bistable potential (also called stochastic Cahn- 
Allen or model A) are analyzed for initial data that are close to profiles that connect the 
two phases. It is shownthat the location of the phase boundary performs a Brownian 
motion. These results have been improved in a number of ways, notably to include small 
asymmetries that result in a drift for the arising diffusion process [7j and to deal with 
macroscopically finite volumes [4] (which introduce a repulsive effect approaching the 
boundary). Also the case of stochastic phase field equations has been considered [5]. 

For interacting particle systems results have been obtained for the zero temperature 
limit of d-dimensional Brownian particles interacting via local pair potentials in [16] : in 
this case the frozen clusters perform a Brownian motion and, in one dimension, also the 
merging of clusters is analyzed [TTJ. In this case the very small temperature is the small 
noise from which cluster diffusion originates. With respect to |16[ I17j. our results hold for 
any super-critical interaction, but of course our system is of mean field type. It is also 
interesting to observe that for the model in |16[ [T7] establishing the stability of the frozen 
clusters is the crucial issue, because the motion of the center of mass is a martingale, i.e. 
there is no drift. A substantial part of our work is in controlling that the drift of the 
center of synchronization vanishes (and controlling the drift is a substantial part also of 
13 [3 [H [5]). This is directly related to the content of § O 

As a matter of fact, in spite of the fact that our work deals directly with an interacting 
system, and not with an SPDE model, our approach is closer to the one in the SPDE 
literature. However, as we have already pointed out, a non negligible point is that we are 
forced to perform an analysis in distribution spaces, in fact Sobolev spaces with negative 
exponent, in contrast to the approach in the space of continuous functions in [8| [7[ |4"[ [5] . 
We point out that approaches to dynamical mean field type systems via Hilbert spaces of 
distribution has been already taken up in [H] but in our case the specific use of weighted 
Sobolev spaces is not only a technical tool, but it is intimately related to the geometry of 
the contractive invariant manifold M. In this sense and because of the iterative procedure 
we apply - originally introduced in [8] - our work is a natural development of [HI 0] . 

An important issue about our model that we have not stressed at all is that propagation 
of chaos holds (see e.g. p2]), in the sense that if the initial condition is given by a product 
measure, then this property is approximately preserved, at least for finite times. Recently 
much work has been done toward establishing quantitative estimates of chaos propagation 
(see for example the references in [H]). On the other hand, like for the model in jllj, we 
know that, for our model, chaos propagation eventually breaks down: this is just because 
one can show by Large Deviations arguments that the empirical measure at equilibrium 
converges in law as N — > oo to the random probability density qx{~), with X a uniform 
random variable on S. But using Theorem 11.11 one can go much farther and show that 
chaos propagation breaks down at times proportional to N. From Theorem 11.11 one can 
actually extract also an accurate description of how the correlations build up due to the 
random motion on M. 
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It is natural to ask whether the type of results we have proven extend to the case in which 
random natural frequencies are present, that is to the disordered version of the model we 
consider that goes under the name of Kuramoto model. The question is natural because 
for the limit PDE [131 [26] there is a contractive manifold similar to M [20] . However the 
results in [27] suggest that a nontrivial dynamics on the contractive manifold is observed 
rather on times proportional to yN and one expects a dynamics with a nontrivial random 
drift. The role of disorder in this type of models is not fully elucidated (see however [12] 
on the critical case) and the global long time dynamics represents a challenging issue. 

The paper is organized as follows: we start off (Section [2]) by introducing the precise 
mathematical set-up and a number of technical results. This will allow us to present 
quantitative heuristic arguments and sketch of poofs. In Section [3] we prove that if the 
system is close to M, it stays so for a long time. We then move on to analyzing the 
dynamics on M (Section^]) and it is here that we show that the drift is negligible. Section[5] 
provides the estimates that guarantee that we do approach M and in Section [6] we collect 
all these estimates and complete the proof our main result (Theorem II. ip . 

2. More on the mathematical set-up and sketch of proofs 

2.1. On the linearized evolution. We introduce the Hilbert space H_i\i q or, more 
generally, the space H-\ jW for a general weight w G C 1 (S;(0,oo)) by using the rigged 
Hilbert space structure [10] with pivot space Lq := {u G L 2 : J s u = 0}. In this way given 
an Hilbert space V C Lq, V dense in Lq, for which the canonical injection of V into Lq is 
continuous, one automatically obtains a representation of V 1 - the dual space - in terms 
of a third Hilbert space into which Lq is canonically and densely injected. If V is the 
closure of {u G C 1 (S;R) : J u = 0} under the squared norm f s (u') 2 /w, that is Hi i/ W , the 
third Hilbert space is precisely JJ_i jt0 . The duality between H\^/ w and H_\ w is denoted 
in principle by ( • , • )h x 1/w ,H-i >w , but less cumbersome notations will be introduced when 
the duality is needed (for example, below we drop the subscripts). 
It is not difficult to see that for u, v G H-x >w 

(u,v)_ hw = fwUV, (2.1) 

where U, respectively V, is the primitive of u (resp. v) such that J s wU = (resp. 
/ § wV = 0), see [6J § 2.2]. More precisely, u G fl-i,™ if there exists U G L 2 (S;R) such 
that JgUvu = and (u,h) = — J^Uh! for every h G Hi ij w . One sees directly also that by 
changing w one produces equivalent H\ w norms [22| §2.1] so, when the geometry of the 
Hilbert space is not crucial, one can simply replace the weight by 1, and in this case we 
simply write H—\. Occasionally we will need also -ff_2 which is introduced in an absolutely 
analogous way. 

Remark 2.1. One observation that is of help in estimating weighted H-i norms is that 
computing the norm of u requires access to IA: in practice if one identifies a primitive 
U of u, then \\u\\ 2 _ lw < LlA 2 w. This is just because U = U + c for some c G R and 
f § U 2 w = f s U 2 w + c 2 f § w. 

The reason for introducing weighted H_i spaces is because, as one can readily verify, 
L q , given in (jl.8p . is symmetric in H_i i/ q . A deeper analysis (cf. [6]) shows that L q is 
essentially self-adjoint, with compact resolvent. The spectrum of — L q lies in [0,oo), there 



RANDOM LONG TIME DYNAMICS FOR MEAN FIELD PLANE ROTATORS 



9 



is an eigenvalue Ao = with one dimensional eigenspace generated by q' . We therefore 
denote the set of eigenvalues of — L q as {Ao,Ai,...}, with Ai > and Xj+i > Xj for 
j = 1,2,.... The set of eigenfunctions is denoted by {&j}j=o,i,... and let us point out 
that it is straightforward to see that ej G C°°(S;M). Moreover, if u G C 2 (8;R) is even 
(respectively, odd), then L q u is even (respectively, odd): the notion of parity is of course 
the one obtained by observing that u G C 2 (S;M) can be extended to a periodic function 
in C 2 (M;R). This implies that one can choose {ej}j=o,i,... with ej that is either even or 
odd, and we will do so. 

Remark 2.2. By rotation symmetry the eigenvalues do not depend on the choice of q(-) G 
M , but the eigenfunctions do depend on it, even if in a rather trivial way: the eigenfunction 
of L q ^ and L q ^, just differ by a rotation ofip' — ip. We will often need to be precise about 
the choice of q(-) and for this it is worthwhile to introduce the notations 

:= L qi> and — L^e^j = Xje^j . (2.2) 

The eigenfunctions are normalized in H_i^/ q ^. 



Remark 2.3. Some expressions involving weighted H_\ norms can be worked out ex- 
plicitly. For example a recurrent expression in what follows is (u,q') ll / q , for u £ H-\ 
and q G M. IfU is the primitive of u such that J s U/q = 0, then we have (u, g')i,i/g = 
j s U(q — c)/q = J^U, where c is uniquely defined by L(q — c)/q = 0, but of course the 
explicit value of c is not used in the final expression. In practice however it may be more 
straightforward to use an arbitrary primitive IA of u (i.e. j^U/q is not necessarily zero) 
for which we have 

(m')i,i/„ = (u(l--\ . (2.3) 



(2.4) 



'1,1/9 

Since now c appears, let us make it explicit: 

2vr 1 



f s l/q 2-kII(2Ki 



2.2. About the manifold M. As we have anticipated, we look at the set of stationary 
solutions M, defined in (|2.2p . as a manifold. For this we introduce 

H-i := L : p--L G #_il , (2.5) 



2vr 

which is a metric space equipped with the distance inherited from -£/_i, that is dist(/zi, /^) = 
Wl^i — 1*2 1| -l- We have M C H_\ and M can be viewed as a smooth one dimensional man- 
ifold in H—\. The tangent space at q G M is q'M and for every u G H-i we define the 
projection P° on this tangent space as P°u = (u, q')-i^/ q q' /(q' , q')-i,i/ q - The following 
result is proven in [351 p. 501] (see also Lemma 5.1]): 

Lemma 2.4. There exists a > such that for all p G N a with 

N a := U qeAI G H-t : ||/x - g||_i < a} , (2.6) 

there is one and only one q =: v(fi) G M such that (/i — 9,90-1,1/9 = 0- Furthermore, the 
mapping fi h-» v(n) is in C°°(iT_i, -ff-i), and (with D the Frechet derivative) 

Dv(v) = I»m ■ (2.7) 
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Note that the empirical (probability) measure fiN,t that describes our system at time t 
is in H_\ (see Appendix [A} and Lemma 12.41 guarantees in particular that as soon as it is 
sufficiently close to M there is a well defined projection v (nN,t) on the manifold. Since the 
manifold is isomorphic to S it is practical to introduce, for /i G H—i, also p(/u) 6 §, uniquely 
defined by v(fj.) = Q , p (^)- It is immediate to see that the projection p is C°°(ll_i, §). 

2.3. A quantitative heuristic analysis: the diffusion coefficient. The proof of The- 
orem [1J] is naturally split into two parts: the approach to M and the motion on M. The 
approach to M is based on the properties of the PDE ([1.40 : in |21j it is shown, using the 
gradient flow structure of ([1.40 , that if the initial condition is not on the stable manifold 
U (see (|1.7|) ) of the unstable stationary solution then the solution converges for time 
going to infinity to one of the probability densities q = q^ S M (of course ip is a function 
of the initial condition), so given a neighborhood of q^ after a finite time (how large it 
depends only on the initial condition), it gets to the chosen neighborhood: due to the 
regularizing properties of the PDE, such a neighborhood can be even in a topology that 
controls all the derivatives |21j . but here there is no point to use a strong topology, since at 
the level of interacting diffusions we deal with a measure (that we inject into i?_i). And 
in fact we have to estimate the distance between the empirical measure and the solution 
to ([1.40 - controlling thus the effect of the noise - but this type of estimates on finite time 
intervals is standard. However here there is a subtle point: the result we are after is a 
matter of fluctuations and it will not come as a surprise that the empirical measure ap- 
proaches M but does not reach it (of course: M just contains smooth functions, and HN,t 
is not a function), but it will stay in a A~ ^-neighborhood (measured in the H-i norm). 
How long will it take to reach such a neighborhood? The approach to M is actually expo- 
nential and driven by the spectral gap (Ai) of the linearized evolution operator (at least 
close to M). Therefore in order to enter such a 7V _1 / 2 -neighborhood a time proportional 
to log A'" appears to be needed, as the quick observation that exp(— \\t) = 0{N~ 1 / 2 ) for 
t > log A/(2Ai) suggests. The proofs on this stage of the evolution are in Section [5j here 
we just stress that 

(1) controlling the effect of the noise on the system on times 0(log N) is in any case 
sensibly easier than controlling it on times of order N, which is our final aim; 

(2) on times of order N it is no longer a matter of showing that the empirical measure 
stays close to the solution of the PDE: on such a time scale the noise takes over 
and the finite N system, which has a non-trivial (random) dynamics, substantially 
deviates from the behavior of the solution to the PDE, which just converges to one 
of the stationary profiles. 

Let us therefore assume that the empirical measure is in a iV~ ^-neighborhood of a 
given q = qy,. It is reasonable to assume that the dominating part of the dynamics close to 
q is captured by the operator L q and we want to understand the action of the semigroup 
generated by L q on the noise that stirs the system, on long times. Note that we cannot 
choose arbitrarily long times, in particular not times proportional to N right away, because 
in view of the result we are after the stationary profile q around which we linearize changes 
of an order one amount. We will actually choose some intermediate time scale TV 1 / 10 as 
we will see in § 12.41 and Remark \2.6\ that guarantees that working with L q makes sense, 
i.e. that the projection of the empirical measure on M is still sufficiently close to q. The 
point is that the effect of the noise on intermediate times is very different in the tangential 
direction and the orthogonal directions to M, simply because in the orthogonal direction 
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there is a damping, that is absent in the tangential direction. So on intermediate times the 
the leading term in the evolution of the empirical measure turns out to be the projection 
of the evolution on the tangential direction, that is (q' ,fJ,N,t — <2 , )-i,i/g/||<Z / ||-i,i/g- O ne can 
now use Remark 12,31 to obtain 



{q, m,t - q)_ 1>1/q = -J K.{0) (» N , t ( d9) - q{9) AO) , (2.8) 

with K, a primitive of 1— c/q (c given in Remark l2.3p . By applying Ito's formula we see that 
the term in (|2.8p can be written as the sum of a drift term and of a martingale term. It is 
not difficult to see that to leading order the drift term is zero (a more attentive analysis 
shows that one has to show that the next order correction does not give a contribution, 
but we come back to this below). The quadratic variation of the martingale term instead 
turns out to be equal to t/N times 



(fC\e)) 2 q(9) M = l- S^- = \\q'\\\ 1/q • (2.9) 

Js i /q 



Since qy>+ £ = Qip — sq'^ + • • • (note that q'^ is not normalized) , (12. 9ft suggests that the 
diffusion coefficient Dk in Therem ll.il is ||</||Ii x / q , which coincides with (jl.lip . 

To make this procedure work one has to carefully put together the analysis on the 
intermediate time scale, by setting up an adequate iterative scheme. Several delicate 
issues arise and one of the challenging points is precisely to control that the drift can be 
neglected. In fact the first order expansion of the projection that we have used 

p (g* + h) = * - , (M ?~ 1 ' 1/g + o(IHI 2 _i) , (2.10) 
[q > q )-i,i/q 

is not accurate enough and one has to go to the next order, see Lemma [A,5[ This is due to 
the fact that the random contribution, which in principle appears as first order, fluctuates 
and generates a cancellation, so in the end the term is of second order. 



Remark 2.5. It is natural to expect that Theorem M . 1\ holds true also whenpo G U and this 
is just because the evolution is attracted to ^- and then the noise will cause an escape from 
this unstable profile after a time oc logiV, since the exponential instability will make the 
fluctuations grow exponentially with a rate which is just given by the linearized dynamics 
(linearized around ^- of course). Arguments in this spirit can be found for example in \30\ 
Ch. 5], see [3] and references therein for the finite dimensional counterpart. However 

(1) this is not so straightforward because it requires a good control on the dynamics on 
and around the heteroclinic orbits linking ^- to Mq \21\ Section 5]; 

(2) the statement would require more details about the initial condition: the simple 
convergence to a point on U is largely non sufficient ( the fluctuations of the initial 
conditions now matter!); 

(3) in general the initial phase ipo on M is certainly going to be random: if the ini- 
tial condition is rotation invariant (at least in law), like if {(pjj }f=i ... ,jv « r e IID 
variables uniformly distributed on S or if ip^ N = 2-kj/N , one expects tpo to be uni- 
formly distributed on S. Note however that uniform distribution of ipo is definitely 
not expected in the general case and asymmetries in the initial condition should 
affect the distribution ofipQ. 
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2.4. The iterative scheme. As we have explained in § 12.31 t ne analysis close to M 
requires an iterative procedure, which we introduce here. We assume that at t = the 
system is already close to M, while in practice this will happen after some time: in 
Section [6] we explain how to put together the results on the early stage of the evolution 
and the analysis close to M, that we start here. So, for = /i7v,o = jr Ylf=i $ ^ N such 
that dist(/xo, M) < a (here and below dist(-, •) is the distance built with the norm of 
by Lemma EH we can define ipo = p(/^o)- Applying the ltd formula to v t = [i t — q^ , we 
see that 



[ e- {t - s ^d e [v s J*v s \ds + Z t , (2.11) 
Jo 



where 

Z < = n E / d °' Q t-s (0, <d' N ) dWi , (2.12) 

3=1 J ° 

and Qf° '(6,9') is the kernel of e~ sL ^o in L 2 . The evolution equation (12. Ill) and the noise 
term (12. 121) have a meaning in H_i, as well as the recentered empirical measures Uf, and it 
is in this sense that we will use them: we detail this in Appendix [Aj where one finds also 
an explicit expression and some basic facts about the kernel Qf°(9,9'). We have started 
here an abuse of notation that will be persistent through the text: de>Qf° s \ 9, fi' N ^j stands 

totdpgt. (MO 

Equations (|2. 1 1 p — (|2. 12|) are useful tools as long as we can properly define the phase 
associated to the empirical measure of the system and that this phase is close to ipo: in 
view of the result we want to prove, this is expected to be true for a long time, but it is 
certainly expected to fail for times of the order of N, since on this timescale the phase 
does change of an amount that does not vanish as TV" becomes large. 

The idea is therefore to divide the evolution of the particle system up to a final time 

proportional to iV into n = njy N -^^ oo time intervals [Tj,Tj + i], where Tj = iT and 
T = T(N) is chosen close to a fractional power of iV (see Remark 12. 6p . Moreover i runs 
from 1 up to n = tin so that tinT(N) = T nN and liniTv T nN /N is equal to a positive 
constant (the rj of Theorem II. lj) . If the empirical measure jit stays close to the manifold 
M, we can define the projections of fj,T h and successively update the reentering phase at 
all times T^. The point then will be essentially to show that the process given by these 
phases, on the time scale oc N, converges to a Brownian motion. 

More formally, we construct the following iterative scheme: we choose 

a = a N :=\N 2 < y / T/N] N ^0, (2.13) 

( > (see Remark 12. 6p . we set r° = and for k = 1, 2, ... we define 

V-fc-i := pOuiw), (2.14) 

if dist()UT fe _i) M) < &n and 

t «n = 4 ll {vV<T )t - 1 } + inf ^ G [Tk-i,T k ], H^-g^.J-i > a N }l {r k- 1>Tk i} . (2.15) 
Then we set 

^ := fit - , (2.16) 

for t E [Tfc_i,Tfc] and t < t*L, and otherwise := v \ for every t > t£ (of course 
can be smaller than T^-i and, in this case, the definition becomes redundant). Therefore 
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the v process we have just defined solves for t G \Tk-i-, Tk 



-(tA7 



rtAr 



where 



N 

V 



.7=1 • /Jfe - 1 



(2.17) 
(2.18) 



Once again, we refer to Appendix [A] for the precise meaning of ()2. 17[) and (|2.18p . 

Remark 2.6. For the remainder of the paper we choose T(N) ~ N 1 / 10 and ( < 1/100. 

The two exponents do not have any particular meaning: a look at the argument shows that 
the exponent for T(N) has in any case to be chosen smaller than 1/2, but then a number 
of technical estimates enter the game and we have settled for a value 1/10 without trying 
to get the optimal value that comes out of the method we use. 

3. A PRIORI ESTIMATES: PERSISTENCE OF PROXIMITY TO M 

The aim of this section is to prove that, if we are (say, at time zero) sufficiently close to 
M, we stay close to M for times O(N). The arguments in this section justify the choice 
of the proximity parameter o~n that we have made in the iterative scheme. We first prove 
some estimates on the size of the noise term and then we will give the estimates on the 
empirical measure. 



3.1. Noise estimates. We define the event 



B 



N 



sup sup Z t 



1-1 



n 



sup sup 

HK»-l(6[T k ,T H1 ] 



Z 



-i 



AK , (3.1) 



where Z k,J ~ is defined precisely like Z^, see (|2.18|) . except for the replacement of Gf_ s (-, ■ 



with gf_ s (9,0') - ^_ 1 ,o(«)/ 4 _ 1 ,o(«| 



Lemma 3.1. lim 



(B 



1. 



Proof. In order to perform the estimates we introduce and work with approximated ver- 
sions of Z t fc and Z k ' x (see Lemma |A.4[) . Define for T^—i < t' < t 

Z lt> = M E / d ^ts X (*> dWi . (3.2) 

The kernel Q^' 1 in this case is (cf. Appendix [A]) 

oo 
1=0 
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where A/ are the ordered eigenvalues of — L^ k _ 1 , e^ k _ 1: i are the associated eigenfunctions 
of unit norm in H^i/q^ , cf. Remark 12.21 and f\p k _ 1 ,i are the eigenfunctions of ^ 

the adjoint in L 2 (see Appendix |A|) , 

Very much in in the same way we define 



yk,A- 



4E/ d el gt--\e^ N )dwi 



with 



(3.4) 



(3.5) 



i=i 



We decompose for Tk-i < s' < s < t and s' < t' < t 

N 



ryk iyk 



j=l ^^k—i 



+ ^Y.j s , d ^t-u(0^ N )dW^ (3.6) 



and an absolutely analogous formula holds for Z k,± : in fact the bounds for Z k and Z k,± 
are obtained with the same technique even if the results are slightly different due to the 
presence of the zero eigenvalue in Z k . Moreover we apply ||a + 6|| 2 ^ 2(||a|| 2 + ||6|| 2 ),so 
that we can estimate the two terms in the right-hand side of (13. 6ft separately. 

And we start with the second term of the right-hand side in (|3.6p : by the orthogonality 
properties of the eigenvectors we obtain 



N rt 



^E/ d e ,g^(,rt N )awt 

3 = 1 Js ' 



-1,1/9 



. oo N rt > r t' 

= E E 

1=0 j,j'=l 

and by taking the expectation 



l f^M )fk-uiW ) dW l dW l> . ( 3 - 7 ) 



E 



N „+< 



7 = 1 Js ' 



-1,1/9. 



1 

iV 2 



oo N r t 



EE 

1=0 j=l 

By Corollary IB. 61 there exists a constant C\ such that 

2 



2(t-u)X, 



E 



(fk->M N )) 2 ] du - ( 3 - 8 ) 



E 



i N ft' v ft' 

Jr E I, d ^tu l (•> rt N ) dw* ^ ^ E /, e " 2(t "" )Ai d - • ( 3 - 9 ) 



-1,1/3. 



z=o 
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Proposition IB.41 Remark IB, 3\ leads us to 

oo ft ' 00 rt' ,2 00 1 / 

J- / e -2(t-u)x ldu ^ e-MW« < (1 



(=0 



1=0 



1=0 



-{t'-s% 



(3.10) 



where the addend with I = (times C) has to be read as t' — s' . The right-most term in 
rf3~T0]) for t' - s' > 1 can be bounded by C(t' - s') + CJ2Zi V* 2 < 3C7(tf - s'). Instead 
for f - s' < 1 we decompose the same term and then estimate as follows: 



L_(tf- a ')-V2j 



1 



z=o 



/= L (t'-s') _1/2 -l+l 
L(t'- S ')" 1/2 J OO 

< E (*'-*') + E 



-(t'- s ')t 



z=o 



«=L(t'-s')" 1/2 J+l 



-2 < (3 + 2C)v / ^7, (3.11) 



where for the first term we have used (1 — exp(— a)) < a, for a > 0. Therefore we 
have proven that there exists C such that for every k and every s, s', t, t' such that 
T/j_i < s' < s < t and s' < t' < t we have 



E 



^E /, ^S& l (-,^)d^ 

7 = 1 ^ 



-1,1/9 



N 



with 



hi(u) : = u 1/2 l [0)1 )(n) + ul^ 



00) 



(3.12) 



(3.13) 



We can do better in the case of Q^- 1 ^, for which a direct inspection of the argument 
we have just presented shows that the linearly growing term in the estimate can be avoided 
(since the term I = is no longer there) and the net result is 



E 



1 N r* 

^E / ft'Sfc 1 ' (-^u N )dwi 

j= i J • 



-1,1/9. 



<: 



Ch 2 (t> - s') 
N 



with ho defined as 



h 2 (u) = n 1/2 l [0il) (u) + l [lj0o) . 
For what concerns the first term in the right-hand side of (|3.6[) . we have 



E 



1 N /■«' / 

E j T (d^i-u (■. ^) - (•> ^)) 



-1,1/9. 



00 N 



1=1 j=l " T k-i 



MEEL (e-<^-e-<— E (4_ 1 ,(^ 7V )) 2 d«, (3.16) 



(3.14) 



(3.15) 
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and, by proceeding like for (|3.9|) . we see that the expression in (|3.16p is bounded by 

2 



Ci 
N 



oo 1 



1=1 



-kits') 



A, 



<ts')b - 1 



N 



i=i 



I 2 



(3.17) 



This last term is estimated once again by separating the two cases of t — s' small and large. 
The net result is that there exists C > such that for every k, every s, s' and t such that 
Tfc < s' < s < t we have 



E 



1 - 

-Y 



N 



3=1 



de,GtX\-^)-d e .Gtu 1 (-^ N )) 



-1,1/9. 



^ Ch 2 (t-s'). 



(3.18) 



In order to complete the proof of Lemma [3.11 quadratic estimates do not suffice: we need 
to generalize (|3.12p . (13.14p and (13.18j) to larger exponents. We actually need estimates on 
moments of order 2m, with m finite, but sufficiently large, so to apply the standard 
Kolmogorov Lemma type estimates and get uniform bounds. We are going to use 



|o + 6|| m < m(\\a\ 



(3.19) 



but actually we will not track the m dependence of the constants. We aim at showing 
that the expectation of the moments of order 2m of the quantities we are interested in 
are bounded by the m th power of the estimate we found in the quadratic case, times an 
m-dependent constant. 

For m = 1, 2, . . ., the m th -power of the expression in (|3.7p gives 



1 N rt' 
j i J ' 



2m 



-1,1/9 



1 



N 



h,... ,l m =0j 1 ,j' v ...,j m ,j' m =l 

in which we have introduced the random variables 

r-t' 

p -\i(t-u) f/ 



(3.20) 



(3.21) 



We now take the expectation of both terms in ()3.20p and all the terms in the sum that 
do not include an even number of each Brownian motion vanish. The number of non-zero 
terms in the expectation can thus be bounded by (2m)\N m . Applying the Ltd formula 
to each of these non-zero terms, we get at most (2m)!/(2 m m!) terms (the number of 
possibilities classifying 2m elements in couples) of the type I\- ■ ■ I m , where 



h — h(hih) - 
We now observe that 



du . 



(3.22) 
(3.23) 
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and for each index index in the first sum in the right-hand side of (|3.2Up gives rise 
either directly to a term Ik(h) (for this it is needed that the two terms share the Brownian 
motion), or in the arising products the terms Ik(k,h') are associated with a term of 
the type Ik(h,h")- Therefore the expression obtained after applying ltd formula can be 
bounded by a sum of terms of the type l\ • • • I m , with 

-t' 



-2Ai(t-u)]g 



du . 



(3.24) 



Therefore we are facing the same estimates that we have encountered in the quadratic 
case, see (|3.8p and (|3,9p . except of course for combinatorial contribution. In the end we 
obtain that there exists C = C m such that 

2m 



E 
E 

In a similar way 

N 



£ 

3=1 

N ,f 



< c 



hf{i! - s') 



^E I de'9tu'\^ N )^l 



-1,1/9. 
2m 

-1,1/9 



^ C 



2m 



N r - 



hf(t' 



NT, 



(3.25) 



(3.26) 



A 1 JTu_ 1 



i=i u±k ~ 1 



-1,1/9 



N 2 ' 



1 oo N 

^ E E G£(s,t,s>)G%(s,t,s>)...Gt(s,t,s>)G%(s,t,s>) (3.27) 



h ,...,lm=0 jl J'l ,—,jm Jm=l 



with 



Gi(s,t,s') 



(3.28) 



We reduce the problem as above to the study of products of integral terms J\ - ■ ■ Jk with 



Ju 



and then, like before, in terms of products of diagonal terms of the type 

2 



Jk 



du . 



(3.29) 



(3.30) 



Again we are reduced to the estimating terms that have already appeared in the quadratic 
case, see (|3.16|) . so we obtain that there exists C = C m such that 



E 



1 N r s ' / n 



2ni 



-1,1/9. 



iC^j-'K (3.31) 
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We now let t' /* t and s' /* s and by applying Fatou's Lemma and Lemma \AA\ 
from([3SD, (13^51) . (13T26D and (|33Tjl we get 

,K{t-s) 



E 







2m" 




z k t-z k s 


-1 







and 



E 



ryk 



2m 

-1 



^ c- 



^ c 



]\fm ' 

h?(t - s) 
N m 



(3.32) 
(3.33) 



The fact that we are allowed to drop the weight in the H_\ norm is of course due to the 
norm equivalence. 

We are now in good condition to apply the Garsia-Rodemich-Rumsey Lemma |33j : 



Lemma 3.2. Let p and \& be continuous, strictly increasing functions on (0, oo) such that 
p(0) = ^(O) = and limt/oo \P(f) = oo. Given T > and 4> continuous on (0, T) and 
taking its values in a Banach space {E, \\.\\), if 

r-T f-T 



J 



p(|t-a|) 



dsdt < B < oo, 



i/ien /or ^ s ^ t ^ T: 



t-s 



1 w 



p( dtt) . 



(3.34) 



(3.35) 



it 2m and ^(u) 



,2m 



(3.36) 



We apply Lemma 13.21 with 

Ht) = zt Tk ^ , P («) 

and Q = 1/100 (Remark I2.6[) . With these choices we can find an explicit constant C = 
C(m, £) such that 

|| Z fc _ Z fc||2m < c ^ _ s)CjB ^ (3.37) 

for every s and t such that Tfc_i ^ s < t ^ and -B is a positive random variable such 
that 



E[B] < 



C 
7V m 



JO 



h?(\t-s\] 
\t-s\ 2 +^ 



ds dt 



(3.38) 



where C is the constant in (|3.32p . For m > 4 the function 1 1— >■ h™(t)/t 2+ ^ , defined for i > 0, 
is increasing (and it tends to zero for i \ 0). So E[£?] is bounded by CN~ m h 7 ^(T)/T^ 
and therefore 



E 



sup 



t ~ ll-l 

It-sK 



rpm—C, 



which leads to 



sup ||Z*||_i ^ a/-A^ 



1 



Then, (recall n = = y) we deduce 



sup sup ||_i ^ 



1 



(3.39) 



(3.40) 



(3.41) 



RANDOM LONG TIME DYNAMICS FOR MEAN FIELD PLANE ROTATORS 



19 



where the right hand side tends to when m is chosen sufficiently large. A similar argument 
gives for 



sup sup \\Z t 

1 < fcsCnTfe.x O^T fe 



-i > -L^ C 



^ c 



NmC-1 • 



(3.42) 
□ 



We now give the main result of the section: 

Proposition 3.3. // ||^o||-i ^ y= and if the event B N defined in (|3.ip is realized (then, 
with probability approaching 1 as N — >• oo) we have 



sup sup 

1 s: k ^nte[T k _ u T k ] 



and 



max 

l < fc n 



-i ^ x/N 



(3.43) 
(3.44) 



Proof. In view of Lemma 13. II we can and will assume that B N is verified. From (I2.17P and 
Lemma I A.2|. we get, for all fc = 1. . . n and t 6 [Tfc_i,Tfc] 

\\v?U < Ce-^^-OH^JI.j + CjT (l + _L=) KY-id S +||^||-i. (3.45) 

The constant C in front of the first term of the right hand side above would be equal to 1 
if we were using the ||.||— 1,1/^, norm. Let us assume that ||z/^ fc Jl-l — N^/y/N, using 
Lemma 13. II we obtain 



_1 Ce- X ^- T ^^ + C(T + VT) sup 



N 

Therefore we readily see that if we define 

t* = sup jt e [T fc _ l5 T fc ] 
we have that for t < t* 



T k -i Ks^t 



Wt\\-1 > 



.fc 1 1 2 



iV 2 < 



+ . (3.46) 



(3.47) 



l^ll-i < C7iV 2C ^ + 2CT 2 iY 4C " 1 + VTN^-^ 



(3.48) 

Therefore since lim^r T 3 A^ _1+4< > = (see Remark l2.6p . for N large enough, we have t* = 
and (|3.43p is reduced to proving n ||z/j, ||— i < N 2C -/\fN for = 1, 2, . . . , n. This holds 
for fc = 1: we are now going to show by induction (|3.44p and therefore that the assumption 
propagates from fc to fc + 1. 

To prove the bound on u^ 1 , assuming the bound on Vq, , we use the smoothness of 

the manifold M. Since we are working in B N , t£ = and we have 



,fc+l 



P, 



P 1 



(3.49) 



P A 
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Since the mapping ip i— > Pi- is smooth on the compact M, we have (cf. 



and the identities 



and 



- tpk-i = p(mtJ - p(m_i) 



(3.50) 

(3.51) 
(3.52) 



combined with the smoothness of p, lead to (using (|3.43p ) 



p-L p-L 



(3.53) 



On the other hand, the smoothness of with respect to ip, (|3.5ip and f|3.52|) imply 



-i 



< c (K-iil-i + Kil-i 



(3.54) 



so the first term in the last line of (|3.49p is of order j^-iV 4 ^, which is much smaller than 
nor N — > oo, since limN TN 2l ^~ 2 = (see Remark 12. 6p . Moreover, Lemma [2.41 implies 



(3.55) 

so the second term in the last line of (|3,49p is also of order jjN^ . Finally, projecting 
(|2.17p on Range ( L q ^ J and by using again Lemma IA.21 we get 



Pi_A\_ x ^ Ce-^- T ^\\4 k J\-i + C (l + ^=) \\v k s\\-lte + 11^11-1, 

(3.56) 

which, since limjv T 4 ^ 1 " 5 ^ = (see Remark EH}, leads for N large enough to 



P- 1 u k 



AT3C/2 



-1 ^ y/N 



(3.57) 



This takes care of the third term in the last line of (|3.49p and by collecting the three 
estimates we obtain (|3.44p and the proof is complete. □ 



4. The effective dynamics on the tangent space 
Proposition 4.1. We have the first order approximation in probability: for every e > 

rk „/ 



£ 

k=l 



(ipk ~ Tpk- 



E 

k=l 



< e 



1. 



(4.1) 
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Proof. Lemma [A. 51 and Proposition 13.31 give (assuming that B N is realized: we will do this 
through all the proof) 



^k ~ ipk- 



90-1,1/9 

(^,(logg^_ 1 ) ,/ )-i,i/ ff ^ , i4 k ^f k _ 1 )-i,i/ q ^ k , (T 

: - + °u 



(4.2) 



2irI$(2Kr) (<?', <?')-i,l/ g 

Since log(^ fe _ 1 ) // is in R(L q ^_ i ), we have 

(*4 fc , (log^^)")-!,!/^., = ((4 k )\ (bg^xn-!,!/^, ( 4 -3) 

and thus using again Proposition 13.31 we get for the second term of the right-hand side 
1 (^(logg^.J")-!,!/^ H.g^)-!.!/^ 



2vr/ 2 (2^r) (<?', <7')-l,l/ 9 



(q',q')-i,i/ q 



-i 



\/iv Viv 

(4-4) 

and hence it is o(T/N), since lirn/v N^/y/T = (see Remark l2.6p . So only the component 
on the tangent space of M at the point ipk-i is of order T/N: 



ipk - V'fc- 



(4.5) 



We now decompose this tangent term. Our goal is to show that the projection of the noise 
Zx k is the only term that gives a non negligible contribution when N goes to infinity. 
However, a direct domination of the remainder - the nonlinear part of the evolution 
equation (|2.17p - using the a priori bound ^ -^=./V 2< > is not sufficient. In fact 



< ^iV 4C • (4.6) 



In order to improve this estimate the strategy is to we re- inject ()2. 17j) into the projection 
(IkiTk)^^^)-!,!^^, where 



(4.7) 



and this leads to a rather long expression 

n n (7, k n' \ , , / n 9 



^2(tpk - i>k-i) = ^2 



k=l 



I (q':Q')-l,l/q 



^2j2 A k,i+o(i), (4. 



k=l i=l 
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with 

-(T k -s)L 



A k) 2 = Is 
^4fc,7 = Is 



iTk - s)L ^ k -ideMs)J*I k (s)] ds,q' y 

ds ><4-! 



(s-T fc _ 1 )L 



z k :j*z k a 



-{T k -s)L q 



-(s-T ft _i)L 



ds '<4-i 
ds, of,, 



-(T fc - S )L g ^_ i 



Z s fe J* e 



-(a-T fc _i)i. 



ds '4=-i 

ds '4-i 



lj5 



Ak ,9 = 1_E 



Z*J*I k {s) ds,^ fc _ 1 , 



(4.9) 



„TAt* 



where we have used the shortcuts E 1 = {t* ^ r^, fe l }, J* + stands for j T " N and and (•, •)* 

The following bound (a direct consequence of Lemma IA.2I and IA.3|) is now going to be 
of help: 



e (Tk s)Lq ^d e [h 1 ( y s)J*h 2 (s)}ds 



^ C 



1 + 



1 



-1 "2 



( S )||-id S . (4.10) 



VTfe - s. 

In fact it is not difficult to see that by using Lemma 13. fl Proposition 13.31 and (I4.10p we 
can efficiently bound all the Akj's, except A k %: 

j>1/2 j.1/2 t^7/2 t^7/2 

|A*, 6 | < ^JV« |A fc>7 | < ^iV 4 <, |A M | < ^iV^ and |A fc , 9 | < ^N<* . 

(4.11) 

Since T^N 9 ^ 1 -> and iV 5 ?/T -)• (see Remark ESD, we get (recall that n = n N = §) 

2^ fc " = 2^ + 2^ A M + o(l) . (4.12) 

fc=i fc=i fc=i 
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For the Af, 3 terms we need to use something more sophisticate. To deal with these 
terms in fact we rely on an averaging phenomena. This method has been used in [1] for 
the same kind of problem. We write the Doob decomposition 



k=l k=l 

where 

Ik = E [Au^Ft^] , 
and M m is a .Fr m -martingale with brackets 

m 

(M) m = Y,( E Ksl^J-Tfc) • 



(4.13) 
(4.14) 
(4.15) 



We have 



Ik = E 



AT2 



dWl 



k=l 



dwi del I 



2vlt{2Kr)q i>k _ 1 {e) 



where dQ$(0,ff) := d e/ gf(6,6'), and from this we obtain 



L T fc ^ Tl, 1 



(4.16) 



Ik = E 



1 



TAr CT , 



ds / ju s (d0') / d0( 1 



2^/ o 2 (2£>)^ fc _ 1 (0) 

e^'a^Ss (0, g')j{9 - o")d e ,gp_-: (e", e') 



with 



1 * 



3=1 



1* ^ Tfc _ x , (4.17) 

<y at K x 



(4.18) 



where {vs'^jsX) is a solution of (jl.ip depending on J r T k _ 1 only through the initial condition 



~j,N 
<Po = 

The stopping time t ctjv is defined as follows: 



T aN := inf{s > 0, \\/j, s 



-1 > a N } . 



(4.19) 



(4.20) 



We now write /Z s (d0') = q$ h _ x (P') d0' + z^(d0') and split the for right hand-side of (liTTTl) 
into the corresponding two terms. 

The term coming <jty fc _j (0') d0' is zero as one can see by using the symmetry: 



gt k - 1 {^k-i + 0,^ + 6') 



(■*/> 



fe-i 



(4.21) 



which follows from the same statement with ipk-i = 0, which in turn is a consequence of 
the representation (jXJ) and of the fact that if ej is even (respectively, odd) then fj is even 
(respectively, odd) too (see Section 12.11 and Section [A} . 
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For the term containing u^(d6') instead we get the bound 



E 



N 



ds / di/f 



kfnl\ 



de / de"{ i 



27rI$(2Kr)q lPk _ 1 (e) 
de'Qp k Zl(0, 0')J{0 - 0")d e ,Gp k -_\(0" , 9') 



< E 



1 

N 



ds\\K 



\H 



\Hi 



(4.22) 



where 
H k s {9') 



s Jn 



de / de"[ i 



2-Kll{2Kr)q^_ l {9) 



do'G^zW, e')j(o - e")d e/ g%zl(e", e 



We now plug in the explicit representation for the kernels: 



(4.23) 



h,h=0 



2tt 



de 



277 



de" i 



2nI*(2Kr)q^_ 1 (e), 
e^ k _ lA (0)J(e - ^)^ fc _ lj / 2 (^)4_ 1 , Jl (^)4_ ll i 2 (^) • (4-24) 



We obtain 

oo 

n^iii < E e 

«,m=0 



-(A ;i +A i2 )(T fc - S ) 



2- 



de 



2tt 



d#" 1 



2^{2Kr)q^_ 1 (e) 



^-i./ilbll/^^i^lloo + ll/^i^llooll/^^lhj • 

(4.25) 

We aim at proving the convergence of this sum. For the integral term, thanks to the 
rotation symetry, we can limit the study to tpk-1 = 0- Since J(6> — e") = —Ksm(e — e") = 
—Ksin{9)cos(e") + K cos{e) sin(#"), we can split these double integrals into products of 
two simple ones. Corollary IB. 51 implies that there exists 1$ in N such that eo,z +2p and 
e (y +2p+i can be writen as 

e<y +2p = pql /2 (ci,i +2 P vi,i +p + c 2 ,i 0+ 2pV2,i +p) + O > ( 4 - 26 ) 

e 0,Z +2p+l = P9o /2 (ci,« +2p+lWl,i +P + c 2,l +2p+lV2,l +p) + > (4.27) 

where 

sup{|c M |,|c 2i i|} < oo (4.28) 

I > h 

and the functions Vn are defined in Proposition IB.4I The Vn are sums and products of 
sines and cosines, and there exists /i£N such that the only non-zero Fourier coefficients 
of Vij 0+P are of index included between h + p — 2 and h + p + 2 and are bounded with 
respect to p. We deduce that the simple integral terms containing eo,z +2p> which are of 
the form 



C 



2- 



±1/2 



(0)e (%(0)d0, 



(4.29) 
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where g is sine or cosine and C is a constant independent of p, are up to a correction 



of order 1/p a bounded linear combination of the Fourier coefficients of q^ 2 or q Q i/z of 
index taken between h + p — 3 and h + p + 3. The same argument applies for eo,/ +2p+i- 
Since these Fourier terms decrease faster than exponentially (this can be seen by observing 
that L exp(acos 9) d6 = 2-7r/ n (a) and that exp(acos(-)) is an entire function), these simple 
integral terms are of order 1/p. Using Remark IB. 31 and Corollary IB.6I we deduce the 
following bound for \\Hg 1 



-1/2 



1- 



\H 



<: C + C y l l+h. e (n-sy^p. + c y e {T k - s )% 
h,h=i 



(4.30) 



i=i 



where the first term of the right hand side corresponds to the case l\ = 0, I2 = in (|4.25p . 
the second term corresponds to l\ > 0, 12 > and the third term to l\ = 0, 12 > or 
h = 0, h > 0. Applying (I4.22P and Proposition 13. 3( we get: 



T i/2 



ds\\H: 



h + i 



1 



+ E12 



■ hi(h 2 + . . 

h,l = l ^ ' /:-:.! 

T 3/2 



< c 



AT3/2 



and thus for iV large enough 

n 
fc=i 

On the other hand, applying Doob Inequality, 



N 2<: , (4.31) 



(4.32) 



sup 

1 < m < n 



\M m \ > y^ 3? 



N 



TN 6 < 



E[(M) n ] < 



and Proposition 

N 1-6C n 

m=l 



it comes 



T 



1,3] 



T 3 



AT1+2C ' 

(4.33) 

Since TN' 1 ^ -)• (see Remark EJjD, the combination of iH2| and (HT331) leads to 

T 



i=l 



0, 



and the proof is complete. 



(4.34) 



□ 



5. Approach to M 

The long time behavior of the solutions to (jl.4p is rather well understood, so, in par- 
ticular, we know that if p$ is not on the set attracted to the unstable solution then it 
converges to one probability density in M (cf. Proposition 15. 2p : this is directly extracted 
from [H[2T]). This takes care of the first stage of the evolution, that is the approach to a 
small neighborhood of M also for the empirical measure. Things however change when the 
empirical measure is at a distance from M which vanishes as — > 00, because the noise 
starts playing a role and the difference between the empirical measure and the solution to 
the Fokker-Planck PDE (jl.4[) is no longer negligible. But we do need to get to distances 
of about A -1 / 2 and this is done by exploiting the approximate contracting properties of 
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the dynamics when the empirical measure is close to qu. We talk about approximate con- 
tracting properties because the noise plays against getting to M and limits the contraction 
effect of the linearized operator. Nevertheless, the proof mimics the deterministic proof of 
nonlinear stability, to which the control of the noise is added. In principle the argument 
is straightforward: one exploits the spectral gap of the linearized evolution. In practice, 
one has to set up an iterative procedure similar to the one developed in Section El because 
the center of synchronization may change somewhat over long times. This procedure is 
however substantially easier than the one presented in Section O mostly because here the 
control required on the noise is for substantially shorter times (logiV versus N\), so we 
will not go through the arguments in full detail again. 

Proposition 5.1. Choose po £ Mi \ U such that (|1 .9|) is satisfied. Then there exists t/'o 
(non random!), that depends on K andpo(-), C, that depends only on K, and a random 
variable ^ n such that 

(.. N 2 <\ 
lim P \\fiN,e N N ~ q* N _i ^ -Jf? I = !> C 5 - 1 ) 

where £jv := [C log N\/N, and Hiun^n = ipo in probability. Moreover for e and en as 
in Theorem \1.1\ we have 

lim P sup ||/ijv,t - QVJLi ^ £ ) =1 i ( 5 - 2 ) 



t£[c N N,e N N] 



Proof. The proof is divided in two parts. First we prove, using the convergence of [it = /J>N,t 
to the deterministic solution pt, that for a given h > (arbitrarily small), there exists to 
such that for e small enough, P(dist(/U7v,t , M) ^ h) — > 1 when N —> oo. Then we show 
that after a time of order logiV, the empirical measure /it moves to a distance ./V^ -1 / 2 
from M, without a macroscopic change of the phase. 

The first part of the proof relies on the following result: 

Proposition 5.2. If po £ Ai\ \ U then there exists if) £ § such that lim^oo Pt = Qip in 
C k (S;R) (for every k). 

Proposition 15.21 is essentially taken from [21], in the sense that it follows by piecing 
together some results taken from |21j . We give below a proof that of course relies on |21| . 
We point out that the very same result can be proven also by adapting entropy production 
arguments, like in [2]. 

Proposition 15. 21 guarantees that the deterministic solution p t converges to a element q^ 
of M. Therefore for t > to, we have that pt is no farther than h/2 from q^ (this is a 
statement that can be made for example in C k , but here we just need it in H^\). Actually, 
it is not difficult to see that one can choose to = ~j^^°Sh, for h sufficiently small (Ai is 
the spectral gat of L q ^ o ), but this is of little relevance here. Applying the Ito formula 



Ht-Pt = e* 2 (/i - po) - / e (t s) 2 [fi s J * fi s - p s J * p s ] ds + z t 

Jo 



(5.3) 



where 

1 N 
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e* 2 is the semi-group of the Laplacian and T~L is the kernel of e s 2 in L 2 . Define Wn = 
{w, Who — Poll-l ^ £}• Using the classical estimate ||e* A//2 u||_i ^ -^=||u||_2 and similar 

argument as in Section^ we deduce that there exist events Wn C Wn such that P(Wjv) — > 
1 and that for all outcomes in Wn we have 

sup \\zt\\-i^J%N^. (5.5) 
From now, we restrict ourselves to Wn- From (j5.3[) we get for all t £ [0,to] 

\\fH-pt\\-i < e + C I -^ = \\ lls -p s \\^ 1 ds + \j%N^ . (5.6) 
Jo V* - « V iV 

The Gronwall-Henry inequality (see [32j) implies that there exists 7 > (independent of 
e and iV) such that 

sup|K,- Pt0 |U < ( e + ,/|ArC) e 7io. ( 5. 7) 

t<t \ v iv y 

So for e = ft/4 and N large enough, [|// io — g^ ||_i ^ h on the event Wat. 

To show that we enter a neighborhood of size slightly larger than iV -1 / 2 , it will be 
we set up an iterative scheme. It is very similar to the one given in Section 12.41 
but with times i, bounded with respect to N. This times are chosen such that after each 
iteration, the distance between the empirical measure and M is at least divided by 2. We 
define ho := h and for to ^ 1 

t m := t m -i + — |loga| , (5.8) 
Al 

h m := -ft m -i , (5.9) 

until the index defined by 

m f := inf {to > 1, ft m < N 2 ^ 1 / 2 } _ (5 10 ) 

The constant a above does not depens on iV and will be chosen below. It is now easy to 
check that toj is of order logiV. Then we define tq := to, and for 1 ^to^toj + 1 

^ m _i := p(^i m -i) , (5.H) 

:=M*-i "flfc.-! ( 5 - 12 ) 
if dist(//t m _ 1 , M) ^ a (see Lemma |2~^|) . We consider for 1 ^ to ^ mj the stopping times 

Tm '■= r m-ll{f m _ 1 <t m _ 1 } 

+ inf{sG [i m -i,i m ], ||/i s -^ m _J|_i > a}l {fm _ 1 >tm _ l} . (5.13) 
and the process solution of 

T, m — 1 71 m _i_ 1 ^ 

— 1 {Tm<im.-l} Zy r m + 1 {?m >tm-l} X 

f e-^"-*- 1 ^*-!^ - e- (tA ^- s)i ^- 1 9,[^*^]d S + Z-~ m > ) , (5.14) 
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where 

N 



~ 1 /"* 

Z ™ " E / de'Gt-s' 1 (0, rt N ) dW! . (5.15) 

j — 1 ^ tm — l 

With the same arguments as given in Lemma [3. 1\ we can prove (recall that mj is of order 
log N) that the probability of the event 



n N :=\ snp snp Z? ^ ^ N< \ (5.16) 

lUm^iii/tM-i^KL -1 V iv J 

tends to 1 when iV — )• oo. From now, we assume that Otv is verified. We insist on the 
fact that the generic constants C appearing in the following do not depend on N, and 
if not mentioned do not depend on a. From Lemma IA.2I and (15.14D we get thet for all 
1 ^ m ^ mf, 



\\PTW-i <Ch m - 1 + C(t + ft) sup |CT-i + \ tm ^ N*. (5.17) 

We now prove that for 1 ^ m ^ mj — 1, ||z/£™_ ||_i ^ h m -x implies ^ h m , and 
that ^ h mf -i implies ||z^ /+1 ||-i ^ N 2 ^ 1 / 2 . Define 

s m : = sup{s G [t m -i,t m ], \\u s ||_x ^ /i^-i}- ( 5 - 18 ) 
Then for s < s* m , if 11*4™ _J|-i ^ h m -\ ,we get using (I5.17P 



\\v?U < Cft^ + C(« + v^)^! + yj tm ^ iV C • (5.19) 

Since N 2 ^ 1 / 2 ^ h m _i, we deduce that = t m if /i is small enough. Then using ()5.14|) 
we get 

IKII-i < Cah^ + Ch^ + ^-^ Nt. (5.20) 
Since ^ ah m -i for /iq small enough, it leads us to (recall that h m ^i = 2h m ) 



IKII-i < 4Cah m + yJ tm ^ N*. (5.21) 

Ifm < m/, ^ tm ^- 1 N c - < Ca/i m and thus < 5CWi, m . Ifm = m f , h m < JV^ -1 / 2 

and thus lli^ 71 ||_i ^ SCaiV 2 ^ 1 / 2 . We now have a good control on u+ = q-r , and 

project it with respect to ^ m (writing fit m = q^ to get a bound for 

We use the same decomposition as the proof of Proposition 13.31 

?, m+1 = qr +»r -Q7 = ^7 [«7 + ~ Q7 ] 

= (Vf -Pf )[g 7 + £f-g 7 l + P r L [g 7 - 97 ] + Pf . (5.22) 
Since the projection p is smouth, we get the bound 

$ m -$m-i\ = | P (/0-p(// tmi )| < C-H^-zx^JI-i ^ CK-^JI-i. (5.23) 
But (I5.2ip implies in particular that 

K||_i<C7(l + 4a)/i m _i 1 (5.24) 
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which implies, using also (|5.23|) . 

$m ~ Vwil < 2C(1 + 4a)/i m _! . (5.25) 

Using similar arguments as in the proof of Proposition 13. 31 (using in particular the smouth- 
ness of the projection Pjj~), we see that the two first terms of the right hand side in (15,22j) 
are of order h^—i' More precisely, there exists a constant C'[a] depending in a (increasing 
in a) such that 

^C'[a]h 2 m _ 1 + C\\»ZU- (5-26) 
So, since ||z^||_i ^ 5Cah m for m < rrif and ||-i ^ §Ca.N 2< *~ 1 / 2 , if ho and a are 

small enough we get ^ h m for m < rrif and ||-i ^ N 2i >~ 1 / 2 . 

We have therefore shown that after a time of order log N, the empirical measure comes 
at distance jV 2< > _1 / 2 from from qj . This angle ib m , corresponds to the angle in the 

Proposition [5TTJ So it remains to prove that ip mf converges to tpQ in probability as N goes 
to infinity. We decompose 

mf 

I • (5-27) 

m=l 

We restrict our study on the event fi/v f] Wn, whose probability tends to 1. Since ||^ — 
Qipo II -l ^ h and the projection p is smouth, we get 

|^o - *k\ ^ Ch (5.28) 
and ()5.25p implies (recall ho = h and h m -i = 2h m ) 

$m-^m-i\^C2 1 - m h. (5.29) 

Consequently for C large enough P[|^ m/ — V'ol > Ch] — »jv->oo 0, which completes the proof 
of (15. ip . The bound (15. 2p is much rougher and it follows directly from the argument we 
have used for establishing (15. ip . This completes the proof of Proposition 15.11 □ 

Proof of Proposition \5.2\ The crucial issues are the gradient flow structure of (|1 .4[) and its 
dissipativity properties. The gradient structure of (ll.4p [6] implies that the functional 

F{p) := \ fp(9) logp(0) d9 - — I [ p(0)cos(e-e')p(6')d0d9' , (5.30) 
2 J § 2 J s Jg 

is non increasing along the time evolution. The dissipativity properties proven in [21\ 
Theorem 2.1] show that for every k G N and a > we can find t such that ||pt||c fc < a f° r 
every t > t. Therefore for any k there exists {t n }n=l,2,... such that t n+ \ — t n > 1 and lim n p tn 
exists in C k and we call it p^ . An immediate consequence is that lim n P(pt n ) = P(Poo ) ■ 
But we can go beyond by introducing the semigroup St associated to (|1.4p . by setting 
SfPt = Pt+v- [2U Theorem 2.2] implies the continuity of this semigroup in C k , so that, 
since for t £ [0, 1] we have t n < t n + t < t n+ \, we obtain P(StPoo) = P{Poo)- Therefore 
dtP(StPoo) = 0, but the condition dtP(pt) = 0, for a solution of (jl.4p . directly implies 
that dgPt = 2dg(ptJ * pt), which is the stationarity condition for (jl.4p . Therefore pt is 
either qy,, for some ip, or it coincides with ^ (see (ll.5j) - (ll.6p ). 

Let us point out that if pt n converges to ^ then {pt}t>o itself converges to a-. This is 
just because P (4-) > Piq^), so that if ]iva. n pt> = q^ and Yvm n p tn = i then it suffices to 
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choose n such that T(pv) < J 7 (^r) and m such that t m > t' n to get FiptQ > ^{Pt m ) — 
J 7 (5=), which is impossible. 

So we have seen that either Imit^ooPt = 57 ° r a h limit points are in M. The stronger 
result we need is the convergence also when the limit point is not This result is 
provided by the nonlinear stability result |211 Therem 4.6] which says that if po is in a 
neighborhood of M (the result is proven for a L 2 neighborhood, which is much more than 
what we need here), then there exists tp such that Yvaxt-^ooPt = Qi> i n C k . 

To complete the proof we need to characterize the portion of M\ which is attracted by 
2^: , that is we need to identify the stable manifold of the unstable point with the set U in 
(flTTp . But this is the content of [HJ Proposition 4.4]. □ 



6. Proof of Theorem 11.11 

The proof of Theorem 11.11 relies on the results of the previous sections and on a conver- 
gence argument of the process in the tangent space that we give here. 

Proof of Theorem 11.11 First of all Proposition 15.11 takes care of the evolution up to time 
Nejsr = ClogN and provides an estimate on the closeness of the empirical measure to 
the manifold M that allows to apply directly Proposition 13.31 and then Proposition 14.11 
Note that the iterative scheme that we have set up in Section 12.41 has been presented 
without asking ipo not to be random or not to depend on N. In fact we start the iterative 
scheme at time Nen and from the random phase *&n of Proposition 15.11 that converges 
in probability to the (non random) value ipQ. Of course there is here an abuse of notation 
in the use of tpo, but notice actually that, by the rotation invariance of the system, we 
can actually consider without loss of generality that the empirical measure fJ>N,c\ogN has 
precisely the phase tpQ. Moreover we make a time shift of Nen, so that the phase is ipQ 
at time To = 0. The result in Theorem 11.11 is given for times starting from Nen and not 
Nen, but as stated in Proposition ETJ the empirical measure stays close to q^ in the time 
interval [Nen, Nen]- Therefore we have the finite sequence of times Tq,T%, . . . ,T n , with 
the corresponding phases ipo, ipi, ip n and we define ipt for every t £ [0, T n ] by linear 
interpolation. We assume T n > t/N. 

We then note that, in view of (|3.43p . the control on the phases, see Proposition 14.11 
on the times T\ , T2 , ■ ■ ■ of our iteration scheme suffices not only to control the distance 
between the empirical measure /xjvt an d (M t , in the H—\ norm, for t = but for every 
t £ [0,T n ]. We are now ready to identify the process Wn,- of Theorem ll.il 

' T := — ~D~k — ' ^ ' 

where we recall that r £ [0, T n /N]. We are therefore left with showing that Wn - converges 
to standard Brownian motion. Note that it would be equally possible and maybe more 
natural to define Wn,ti for t > e^r as in the right-hand side of (|6.ip . but with r replaced 
by r — en, and Wn, t = for r G [0,£jv]. In view of the statement we want to prove this 
detail is irrelevant. 

In proving the convergence to Brownian motion we apply Proposition 14.11 and replace 
the process tp. with the cadlag process tpo + Mn,- £ D([0, T n /N];M) defined by 

M N , T := AM N,k, (6.2) 

keN:T k <Nr 
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and 



AM, 



N,k 



( Z T fc; g^._ 1 )-l,l/gy, fc _ 1 
(l',l')-l,l/q 



(6.3) 



It is straightforward to see that Mjv,- is a martingale with respect to the filtration T T := 
F\jT\/t> where T. is the natural filtration of {W 3 N .}j=\ : ... : N'- the martingale is actually 
in L p , for every p, as the moment estimates is Section [3] show. We can now apply the 
Martingale Invariance Principle in the form given by |23|. Corollary 3.24, Ch. VIII] to 
Mat- for continuous time martingales: the hypotheses to verify in the case of piecewise 
constant cadlag martingales boil down to the variance convergence condition that for every 
r G [0,7-/] 



lim 



E 



E 



tD 



K ! 



(AM N , Tk y F Tk _, 

k&N:T k <rN 

in probability, and the Lindeberg condition that for every e > in probability we have 



(6.4) 



lim 



E 



E 



(AM^tJ 2 ; AM 2 NjTk >e T Tk _, 



k£N:T k <rN 

For what concerns (16.41) we have 
E 



0. 



(6.5) 



(AMjv,r fc ) T Th _ x 



4_ 1>O W)%iV lS (d0)d S . (6.6) 



Now take the sum over k and use the uniform estimate (|3.43|) of Proposition 13.31 to re- 
place the empirical measure with q^ T (9)d6. Since a direct computation shows that 

W(0h,o) 2 q46)M = 1, follows" 1 

For what concerns (|6.5p we remark that, by the Markov inequality, it suffices to show 
that 

\4 



lim 



E 



0. 



(6.7) 



E[(AM N , Tk r T Th _ x 

keN:T k <rN 

Actually one can show that there exists a non random constant C such that almost surely 

2 



E 



(AM, 



< C 



(6.8) 



This is an immediate consequence of (|3.32p . but of course, since we are projecting on q' and 
since we are just considering the fourth moment, a similar estimate can be easily obtained 
explicitly by proceeding like for (|6.4p and by using the fact that ||/^ lloo = 1 1 /o I loo < °o- 
Of course (j6.7|) follows from (|6.8p . 

Therefore Mjy r £ D([0, r/];R) converges in law to W./H^H-i where W. is a standard 
Brownian motion. This is almost the result we want (recall that Dk = l/lk'lUi.i/g); since 
we Mn : -/Dk differs from Wat,, just for the fact that they interpolate in a different way 
between the times (where the coincide) and that in the case of Wn . the convergence 
is in C°([0, r/];R). But (j6.8j) guarantees that the sum of the fourth power of the jumps 
of M|\r )t adds up to 0(T 2 /N) = o(l) in probability, so the supremum of the jumps is 
o(l), and therefore the convergence for M^ t . £ D([0, r/];M) implies the convergence of 
Wn,- S C°([0, r/]; R). The proof of Theorem 11.11 is therefore complete. □ 
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Appendix A. The evolution in H_i 

In what follows we fix q in the invariant manifold M (see (II, 6\\ ). Unlike the rest of the 
paper here we do not identify q with q^ and then with ip, so in particular we write L q 

(and not L^), Gf(-) (and not Gf(-) like in (|2.12p ). and so on. We work with the signed 
measure 

v N ,t(M) ■= fi N ,t(d9)-q(e)d9, (A.l) 
which can be seen as an element of H-\. This is simply because it is the difference of 
two probability measures. In fact, if fj, G Mi, 9 i— > fi([0,9]) is a primitive of fi and, by 
Remark EU \\fi - i/|| 2 _ 1 < f s (p([0,6]) - v{[Q,9])) 2 d9 < 2vr. Therefore < v^vf: 

of course this quick argument needs to be cleaned up by first smoothing the measures. 
That is, we introduce an approximate identity (j) n G C°° (4> n > 0, <f> n {6) = for 6 G 
[1/n, 2tt — 1/n], J s n = 1 and lim n J" s F<^ n = F(0) for every F G C°). We then introduce 
the probability density 8 \-t fi n (9) := f s <fr n (9 — 9')fj,(d9') and verify that 

H^n-Mmllii < — -r-j- r, (A. 2) 

mm(n, m) 

so that lim n ^ n exists in H_i (of course the limit exists also weakly and it is fi). 
We aim at proving: 

Proposition A.l. If {<ft' N }t>o j=i ...,N solves (11.11) then . G C°([0, oo); H-\) and we 
have 

v N)t = exp(tL q )u N)0 - / exp((t - s)L q )d ((J * u NiS )u N}S ) ds + Z N)t , (A.3) 

Jo 

where Z]y )t is the limit in H-i as r /~ t of Z^,t,T, where 

1 N r 

Z N ,t, T (9) ■■= TrE/ d e ,Gl s {9,^ N )dWi (A.4) 

3=1 J ° 

Moreover all the terms appearing in the right-hand side of (|A.3p , as functions of time, are 
in C°([0,oo);iT_i). 

Proof For (t,6>) i-> F 4 (6>) in C 1 ' 2 (M+ x S;R), from dTTTJ we directly obtain 

F t (9)v N ,t(d9) = [ F (9)u N)0 (d9) + f [ (L*F S ) [9)v N , s ( dO) ds 
Js Jo Js 

f ! d s F s {9)v N>s (d9)ds+ f f deF 8 (0)(J*v Ni8 )(0)v N , a (M)ds + Z% t , (A.5) 
o Js Jo Js 



+ 

where 



j 

and L* is the adjoint in Lq of L q , that is 



1 r* N 

^ = ^/ e^wU,^, (a.6) 



L > = \v" + (J * q)v' - J * (qv') - [ (J* q)v , (A.7) 
2 Js 



for v G C 2 (§; R) such that / g u = 

We sum up here some useful properties of L*: 
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(1) In [B] it is shown that the LQ-norm is equivalent to the Dirichlet form norm of 
L q : the squared Dirichlet form norm of u is ||n|| 2 1 l , + (u, (—L q )u)_ 1 y . On the 
other hand it is straightforward to see that the properties of L q in fl-^i/g, notably 
the fact that it is self- adjoint and that it has compact resolvent, still hold true in 
the space of the Dirichlet form. So L q has compact resolvent in Lq, which directly 
implies that L* has compact resolvent and the very same spectrum (see e.g. [31\ 
VI.5]). 

(2) Recall that we denote by {ej}j=o,i,... a complete set of eigenvectors of L q which is 
orthonormal in H_iij q and observe that there is a unique solution fj to 

A q f 3 (9) := -dg{q{9)d e fj{9)) = ej (0), (A.8) 

such that Jg fj = 0. More generally, A q is a bijection from {u G C°° : f s u = 0} 
to itself: in fact, v = A q u is equivalent to u' = —V/q in our standard notations, 
which determines u since J s u = 0. In particular /' = —£j/q and fj G C°°, since 
ej is C°°, and one obtains 




(A.9) 



By using the fact that q(-) is even, one verifies directly also that if ej is even 
(respectively, odd) - recall from Section 12.11 that ej is either even or odd - the fj 
is even (respectively, odd) too. 
(3) By observing also that L q A q = A q L* one verifies that {fj}j=o,i,... is a complete set 
of eigenfunctions for L* and, of course, L*fj = —\jfj- 

Therefore for every t > and s < t we can define F s (9) = (exp((t — s)L*)F)(9) for 
F G Lq and standard parabolic regularity [?] results imply that F s (-) is C°° for s < t (in 
our case this can be proven directly by using the Fourier transform, like in |21| . but for 
what follows we choose F G C 2 and the regularity result is even more straightforward). 
By plugging this choice into (jA.5h we obtain 

f F(e)u Ntt (dd) = I \ew(tL* q )F)(0)v N , o (M) 

+ f I d e (exp((t-s)L* q )F)(e)(J*v N>s )(e)u N>s (d9)ds + ZF >t . (A.10) 
Jo Js 

At this point we step to looking at fjsrj as an element of and we reconsider (jA.lOj) 
with this novel viewpoint. 

First of all J s F(9)i'N,t(d9) = {F,UN,t)i,-li where (•, i s the duality between Hi 

and H-i (cf. Sec. 12. ip . For the first term in the right-hand side we observe that, for 
v G H-x w e have (exp(tL*)i ? , = {F,exp(tLq)v)i-\: this is because this relation 

holds when v G Lq (in this case the duality can be replaced by the L 2 scalar product) and 
because one can choose a sequence {f n }n=i,2,...) v n G Lq such that v n — > v in H-\ (one 
can choose v n = 4> n * v) so that 

{exp(tL*)F,v)i i = lim(F,exp(tL q )v n ) 2 = (F,exp(tL q )v)i _i , (A.ll) 

n 

where we have used the continuity properties of the duality and of the semigroup operator. 
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For the second term in the right-hand side of (jA.lOp we write 

I d e {exp((t-s)L* q )F)(e)(J*v NjS )(9)v N>s {d9)ds = 
o Js 

[ ((J * i/jv,s)^exp((t — s)L*)F, ^jv, s )i,-i ds , (A.12) 
Jo 

We now introduce v n ^ s := <f> n * vn,s so that for every s G [0,t) 

((J* VN,s)dexp((t - s)L*)F, vn, s )i ,-i = - lim (F, exp((i - s)L q )d((J * VN,s)v n ,s)) 2 

= -(F,exp((i - s)L q )d((J * vn,s)vn,s))i,-i, 

(A.13) 

where in the last step we have used the fact that exp((t — s)L q ) is a continuous operator 
from H_2 to H_\ (Lemma IA.2|) . Notice moreover that we have 

| (F,exp((t - s)L q )d((J * VN,s)Vn,s)) 2 I < \\4 > n\\l\\d((J*VN,s)dexp((t-s)L q )F)\\ 2 \\vN,s\\-l , 

(A.14) 

and, since J(-) = — -ftTsin(-), one sees that this expression is bounded by a constant times 
|| F" || 2, uniformly in n and s < t. Such a bound tells us that one can exchange limit and 
integration in 

/ lim (F,exp((i - s)L q )d(( J * VN,s)v n ,s)) ds , (A.15) 
Jo n \ J 2 

and then, for fixed n one can of course exchange integral in ds and integral in d9. At 

this point we appeal again to Lemma IA.2I that guarantees that Jq exp((i — s)L q )d((J * 

VN,s)vn,s) ds converges, in to Jjf exp((t — s)L q )d{(J * vn,s)vn,s) ds: note in fact that 

||d((J * zvjv iS )u)||_2 < cj|M|_i so that (by Lemma I A. 21 ) 



< 

-l 



exp((t - s)L q )d((J * v>N,s)(v njS - v n ' iS ) ds 

cjC J ^1 + -j=L=^J \\v ntS - v n ', 8 \\-i ds , (A.16) 
and the right-hand side vanishes for min(n, n') — > oo. Therefore we obtain 



' fl (exp((t - s)L;)F)(9)(J * u N .)(0)v N JM) ds 



J§ 



q 

-t 



(F, [ exp{(t-s)L q )a((J*u N , a )u Nta )ds) lt - 1 . (A.17) 
J o 

We are left with the last term in (|A.10|h It is now useful to use the kernel of the 
Z/ g -semigroup in Lq 

oo 

gi(e,e') ■.= J2 e M-^i)ei(0)fi(0'), (A.is) 

1=0 

so that 

(u,exp(sL q )v) 2 = (exp(sL*)u,v) =11 u{9)G q s (9, 9')v(9') d9d9' . (A.19) 

Js Js 

Note also that, for s > 0, Q q s is C°° in both variables, by the standard parabolic regularity 
results we have mentioned above. So, for every r < t, 9 h-» Zn^, t {8) (recall (IA.4P ) is well 
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defined and smooth in 9. But Lemma lA. 41 tells us that \\m T/ /~ t Zn^,t) exists in H—\. If we 
call the limit Zjv,t) we directly see that (recall (|A.6j) ) 

Z£ <t = (F^nj)!,-!. (A.20) 

Therefore we have shown that (jA.lOp implies the validity of (|A.3p if we take the duality 
with respect to an arbitrary F G C 2 . But we have also shown that every term in (1A.3P is 
in H-i, therefore the equation extends to F £ Hi and (|A.3j) is proven. 

The continuity claimed in the statement follows by the continuity of the three terms 
in the right-hand side of (|A.3p , The continuity of the first term is immediate from the 
properties of the semigroup. The continuity of the second term follows from a direct 
estimate by applying both bounds in Lemma lA.2l Finally the continuity of the third term 
is claimed in Lemma IA.4I The proof of Proposition lA.il is therefore complete. □ 

Lemma A. 2. For r > the operator exp(rL g ) extends to a bounded operator from H_2 
to H-i and there exists C > such that for every r > 

||exp(rL g )u||_ 1 < C (l + -^=\ ||«[|_ a , (A.21) 

and such that for every e G (0, 1/2) we have 

||exp((r + 5)L q )u - exp(rL (? )u||_ 1 < C5 £ (l + -^j^j NI-2 , (A.22) 

for every r > and 5 > 0. 



Proof. We introduce the interpolation spaces associated to L q that is the (Hilbert) spaces 

( oo oo ~| 

V m := \u = ^2u k e k , ]T(l + A fe r*4<ocL (A.23) 



fc=0 k=0 



associated with the norms 



Ul 2 

| t* 1 1 yrn 



||(1 - L q r' 2 u\t = E( X + X ^ u l ■ (A.24) 

k=0 

It is proven in |22[ Remark A.l] that the norms ||.||v n and ||-||n-i are equivalent. This 
equivalence can also be deduced from Remark IB. 3 1 In particular ||-||y-i and ||.||-2 are 
equivalent, so we will prove (|A.2p with ||. and ||.||y-i. For all u = YlT=o u k e k, we 
extend e rLq u as 



e rL "u - ^-- A * 



Ys^Ukek, (A.25) 

fc=0 

and we deduce 

oo 2 

\\e TL «u\\\ 1/q = £(1 + Afc )e- 2A ^-^- . (A.26) 
k=o 1 + Ak 

But if we define f(y) := (1 + y)e~ 2yT , it is easy to see that for all y ^ 0, there exist C 

such that f{y) C 2 (l + -j^j , which with (L\~24|) . and (|A~26l) gives the first inequality. 
For the second inequality we make a similar spectral decomposition and we obtain 



\ e (r+s)L qu _ e rL, u ||2 1 = yVj + Xk )e~ 2X ^ (1 - exp(-,5A fc )) 2 . (A.27) 

. — : 1 + *k 



k=0 
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We then use (1 — exp(-x)) < x e for x > and (1 + x)x 2e exp(— xr) < C 2 (l + r 6 1 / 2 ) 2 , 
for a suitable C which can be chosen independent of e £ (0/1/2). □ 



Lemma A. 3. For all u,v £ there exists C > such that 

\\d e {uJ * v)\\h^ 2 ^ C'IMI-H-iIMIh-i • 



(A.28) 



Proof. For this proof it is practical to write the i^_i-norms by using the Fourier coefficients. 
In fact if u € H_ s , here s = 1 or s = 2, we can define u n = (u,b n ), where b n (6) = 
exp(m#)/27T (note that uq = 0) and 9 *->• X^nez-|n|<Af u " exp(m#) converges as A" — >• oo in 



H s to u. Moreover we have 




1/2 



Since J(0) = — ATsin(^), a direct calculation gives 



de{uJ * v) = Kit 
from which we extract 



(m + l)«i eHm+1) 



\de(uJ * v)\\'_ 



K 2 tt 



^ 4 \m(v-iu m+1 - viu m -i) 



^ Ar 2 7rmax(|w_i| 2 , |t>i| 2 ) 



-2/2 



m [u 



m—l ' "m+l / 



(A.29) 



(A.30) 



< 4AT 2 7r|H| 2 1 ||u|| 2 1 . (A.31) 



□ 



Lemma A. 4. The almost sure limit of Z^ % T as r /* i exists in H_\ and, if we call the 
limit point Z^j, we can choose a continuous version of Z^,-, that is Z^,- £ C°([0, oo); H-\). 

Proof. The claim follows from the same estimates as the one that we have obtained for 
the proof of Lemma 13-H which are however substantially more precise than what we need 
here: recall that now N is fixed, while in Section [3] one of the crucial points is to follow the 
N dependence of the results. Therefore we will not go through the arguments in detail, 
but we just point out that that one goes from Z^,t,r to Z^ t , , see (|3,2p by making obvious 
changes. So, in particular, proceeding like for (|3.25|) we easily gets 

< Ch?(\T-r'\), (A.32) 

where C depends on A~ and m and < r, r' < t. An estimate like ()A.32|) implies almost 
sure Holder continuity of Zfft. , by a direct application of Kolmogorov continuity Lemma 
[33j or by using the Garsia-Rodemich-Rumsey Lemma (Lemma 13. 2p . That is, there exists 
a (positive) random variable X and a positive constant c > such that 

||^JV,t,r — ^JV,t,T'||_ 1 < A|t-t'| c , (A. 33) 

for every < r, r' < t. Therefore the almost sure limit of Zjft T , as r /* t, exists. 



E 



\Zn. 



t,T 



Z 



N,t,T' 



1 2m 
1-1 
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The continuity of the limit follows in the same way, this time using also (|3.31|) . Actually, 
in the proof of Lemma 13,11 we use Lemma IA.4I only to define Z\ as almost sure limit in 
H—\: the proof of continuity is strictly contained in the argument that starts from (|3,32p 
and goes till the end of that proof (but, once again, that proof is substantially more 
informative and involved, since it follows the A-dependence). □ 

A.l. Second order estimates of the projection. As anticipated in $ !2.3l our approach 
requires a control up to and including the second order for the projection map p(-) (recall 
§ 12.21 for the definition). The expansion is with respect to the H-\ distance from the 
manifold M. 

Lemma A. 5. For all q = £ M and h E H—\ with \\h\\—i < a, we have 
{Kq')-x,X/q ( 1 (^,(logg)")-l,l/ g \ 



Proof. For h as in the statement we have that 

{q^ + h- q^ +e , g^ +£ )-i,i/^ +e = , (A.35) 

for e := p{q^ + h) — ijj. Since p{ ) is smooth, we have e = 0(\\h\\_i). By expanding q^ +£ 
with respect to e we see that (fAT35l) impl ICS 

(h + eq'^ - q'^ +e ) = 0(e 3 ) . (A.36) 

Let us rewrite (|A.36P more explicitly (recall Remark 12. 3p as 

L 0> + ■*<»> - i 1 - 2 <(2A»^!w ) d " - o(e3) ■ (A - 37) 

where H is the primitive of h such that J s ^ = 0. At this point we expand also q^+ £ with 
respect to e and, using e = 0(||/i||_i), the parity of q^(- + ip) and Remark 12.31 we get to 

(h, ^)-i,x/ H + g(<, + g 2 7rJ 2 ( 1 2Kr) (1 ° g g)//) ~ 1 - 1 /^ = °(ll /l H-i) • ( A - 38 ) 

Now it suffices to solve this equation for e and perform one last Taylor expansion. □ 

Appendix B. spectral estimates 

The aim of this section is to find approximations of the eigenvalues and eigenfunctions 
of the operators for large eigenvalues. In such a regime we expect the Laplacian to 
dominate and the spectrum of should get close to the one of the Laplacian (as long 
as we deal with large eigenvalues). These are standard estimates, developed for example 
in [28] that we follow, but we could not find in the literature the result for the non-local 
operators we consider. Without loss of generality, we can focus on Lq. We have 

Lqu = -u" — (uJ * qo + qoJ *u)' = -u" — (J *qo)u — (J *%)u — q' J *u — qoJ' *u . (B.l) 

We make a change of variable to get rid of the coefficient of order 1: if we define 

u = = ^y, (B.2) 
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and we observe that ^fq = e J * qo , with J (9) := K cos(9), then we get 

u' = Vooy' + (J* qo)V%y , (B.3) 

u" = y/%y" + 2( J * q )V%y' + (J* q' )V%y + ( J * qofVooy > ( B - 4 ) 
and these two last equations together with (|B,1|) give 

L qo y/qoy, = ^[Vqoy" + 2(J * qo)V%y' + (J * q' )Vqoy + (J * qo) 2 y/qoy] 
- (J * qo)[Vqoy' + (J* qo)Vooy] - (J * q'o)Vqoy - %J * (Vqoy) - qoJ' * (V%y) (b.5) 

which leads, after simplification, to the new operator 

Ly := -zy"-m{y), (B.6) 

where we have set 

m(y) ■= -((J * go) 2 + J * q' )y + -7= J * (Vqoy) + Voo J ' * (yfiov) ■ ( B - 7 ) 



Of course m(y) is a function and when we want to make explicit the ^-dependence we use 
mg(y). Since the operator Lq is negative, we are interested in couples (p,y) solution of 

Ly = -p 2 y , (B.8) 

where p is a positive real number. The method of variation of the parameters shows 
that such solutions exist (for all p > if we do not restrict the study to the 27r-periodic 
eigenfunctions of I) and are of the form 

y (e) = Cl e^ 9 + c 2 e-^ 9 - 4- f ' G(9,e>,p)m e/ (y)de> (B.9) 

V 2/0 Jo 

where 

G(6, 9\ p) = ie^ e - e ') - ie-^0- 6 ') . (B.10) 

We define y± the solution such that c\ = 1, c 2 = 0, and y 2 the solution such that c\ = 0, 
C2 = 1. In what follows we start by getting a first estimate of the eigenfunctions y\ and 
2/2 with respect to p — > 00. This estimate implies a first estimate of the eigenvalue 
—A = — p 2 , and this leads to a new approximation of the eigenfunctions, and thus a new 
approximation of —A. This procedure can be repeated recursively, but for us two steps 
will suffice. 

Lemma B.l. For each A > 0, there exist y\ and y 2 independent (non necessarily periodic) 
eigenfunctions of L associated to —A such that (recall that p=y/\): 

yi (9,p) = e^ + O^j , (B.ll) 

2/ 2 (M = e-^ + oQ , (B.12) 

y[(9,p) = V2pie^ e + 0(1), (B.13) 

y' 2 (e,p) = -V2pie-^ e + 0(l), (B.14) 

where 9 S [0, 2tv] and O(-) is as p tends to infinity (and we stress that here and below the 
O(-) term does not depend on 9 or, equivalently, it is uniform in 9 £ [0,2n]). 
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Proof. We prove the result for y\. The proof for y 2 is similar. We define 



A (9,9',v) = --^-G(e,e',p)mo,(v)l e , <B (B.15) 



so that for 9 e [0,27r], 



m (0) = e ^ i9 + / Ao{d,tf,y x )6ff . (B.16) 
j 

The expression for j/u cf. (|B.9[) . can be iterated arbitrarily many times and it leads to a 
series expression for yi, at least for p sufficiently large. To see this set fo(9) := e^^ 9 and 
observe that 



yi(0 O ) = /o(^o)+V / ••• / A (9 ,9 1 ,A (9 1 ,9 2 ,---A (9 i - 1 ,9 i ,f )---))d9f- d9 m 
j=1 Jo Jo 

+ r f 2W ... r Ao(9o,9 1 ,A (9 1 ,9 2 ,---A {9 m ,9 m+1 ,yi)---))d9 1 --- d9 m+1 . (B.17) 
Jo Jo Jo 

One directly verifies that there exists C = C(K) such that for 9, 9' G [0, 2tt], 

\A (9,9',v)\ < -||«||, (B.18) 

where ||u|| := sup eg [ ,27r] 1^(^)1- Prom (|B. 17[) and using ||/o(-)l = 1 we see that 

^/2^C7\ m /27rC\ m+1 n , ,„ . 

Illftll <! + £(— J + (~J llwlh ( B -!9) 

so for p > 27rC we see that ||yi|| < 00 and we have a series expression for yi, from which 
we directly obtain (|B.llh . 

To deal with y[ we take the derivative of both sides of (|B.9[) with c\ = 1 and c 2 = 0, so 
that 

y[ (9) = v^y^-^T / d e G(9,9',p)m e '(yi)d9'. (B.20) 

V2/0 Jo 

We define the new kernel 

Ax(9,9',v) := --±-d e G(9,9',p)m e '(v)l e , <e , (B.21) 



so we can write 

r 2w 



y[(9) = V2ie^ e + [ * A x (9, 9', yi ) d9' . (B.22) 
J 



Also A\ verifies 

|Ai(0,0',^)| < - sup |u(0)| , (B.23) 

Z 9 6»e[0,27r] 

for a suitable C = C(K) and the same argument as above gives 

-MO) = V2ie^ pie + (-) , (B.24) 
P \PJ 

which is equivalent to (IB. 12ft . □ 
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Lemma B.2. There exists Iq G N such that for all p £ N the eigenvalues of Lq satisfy 

-.2 



A 



Z +2p 



P 



+ 0(vp), 



A/ +2p+l = y + 0(y/p) . 



(B.25) 
(B.26) 



Remark B.3. An immediate consequence of Lemma \B. 2\ and of the basic properties of Lq 
is that there exist C > 1 such that for j = 0, 1, . . . . 

n"2 



J 



(B.27) 



Proof. Let j/i and 7/2 the eigenfunctions of L given by Lemma IB , 1 1 associated to the eigen- 
value —A = — p 2 . As a linear combination of y\ and 2/2 is 2-7r-periodic, the following 
determinant is equal to zero: 

yi(27r)-in(0) 2/2(2^-2/2(0) 
2/1(2^-^(0) 2/2 (2tt) - 2/2(0) 



0. 



Lemma IB. II implies 



e 2V27rpi _ l + Q{ l 
y^piie 2 ^^ - 1) + O(l) 



-2\f2-Kpi _i + Q f I 



-^piie-^P 1 - 1) + 0(1) 



and thus we get 



.rpi _ x|2 = 



We deduce that there exits k £ N such that 



Reciprocally, all p satisfying (|B.3ip satisfies (|B.29|) . so the Lemma follows. 



(B.28) 



(B.29) 



(B.30) 



(B.31) 



□ 



Proposition B.4. There exists I06N such that for allp £ N the eigenvalues of Lq satisfy 

2 ^22 /-JX 
K+2 P = V - + O f - ) , (B.32) 



P 



A 



lo+2p+l 



p- K r A „ ( 1 
+ ( - 



(B.33) 



and any eigenfunction of Lq associated to A; +2p or A/ +2p+i is, up to a correction of order 
1/p 2 , a linear combination of the two functions % vx,i Q +p and V2,i +p, where 

s'm(p6) 



v h i 0+p (9) = cos(p6) 
V2,i +p(®) = sm(p6) + 



p 

cos(p9) 
p 



Kr K 2 r 2 
_ sin (0) + __ sin ( 2 0) 

Kr K 2 r 2 



(B.34) 



From Proposition IB. 41 one can directly extract some important conclusions: let us give 
them before the proof of the proposition. 
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Corollary B.5. There exists Iq G N such that for all p G N and ip G S, £/ie unitary 
(in H_i ±/ q ) eigenf unctions e^ : i 0+ 2 P and e^ j / 0+ 2 P +i of are up to a correction of order 

1/2 

1/p a bounded (with respect to p) linear combination of 9 i— )■ pg, {9)v\^ 0+p {9 — ip) and 



1 /2 

^ P°J {6) v 2,i +p(@ — V") f see Proposition B.J\ for the definition of v\ \ and vi^i). 



Proof. We set tp = without loss of generality. Proposition IB . 41 tells us that the normalized 
eigenfunctions of Lq can be written either as 



sin(pfl) 
cos(p#) 

-,2\ 



if 2 r 2 

_ sin (0) + __ sin ( 2 0) 



+ r p (0)) (B.35) 



where r p (9) = 0(l/p ) and c p is the normalizing constant, or with the analogous expression 
coming from the second line in (|B.34p (but we will deal only with (|B.35P because the other 
case is treated analogously). To estimate c p let us observe that the first two addends in 
(|B.35|) are in H^ 1 (since they are smooth, it suffices to remark that their integral from 
to 2-7T is zero), so r p £ H—i, since the eigenfunction is: of course r p is smooth, since 
the eigenfunction is. Now we claim that the #_i,i/ g norm of cos(p-), that is the first 
addendum, is proportional to 1/p, apart for a correction that is beyond all orders in 1/p, 
while the norm of the two other terms is 0(l/p 2 ). In fact if we set u(9) := cos(p9), then 
U{9) = sin(p9)/p so 

„ „ 1 I f 1- cos(2p9) ,„ _ . 

If we use the standard estimate 



(A;!) 1 / 2 ' 

(B.37) 



Is 



On the other hand 



1 r 1 / \ 2m+k 

'*<*>- Sj[ "< w ^°- Wj, -S rol r(m + * + i) (f) < 

we readily see that 

^l,9 = o('). (B.38) 
q{9) \y/(2pyj 

I "77n d0 = (2vr/ (2Kr)) 2 , (B.39) 

so llitll-i,!/^ * s ec l ua l to c(K)/p, c(K) := \/27rJo(21£Y), up to a correction that decays faster 
than any power of 1/p. 

For the second addendum it suffices to observe that it can be rewritten as a linear 
combination of terms of cos(p'9), with \p — p'\ = 1 and 2. But then the computation 
is very similar to the one that we have done for the first addendum (or, easier, one can 
explicitly compute the H_l norm, without weight). Therefore this term is 0(l/p 2 ). 

For the third addendum we recall that \r p {9)\ < C/p 2 , so that if we set K(0) := 
Jq r p {9') d9' , we have ||7£(#)| < C9/p 2 . Of course 1Z is not necessarily centered, but, by 
using Remark 12 .1\ we see that ||r p ||_i < 2C 2 n 2 /p 2 . 

By collecting the estimates of the three addends we see that 

c p = c{K)p(l + 0{l/p)) , (B.40) 
and this completes the proof of Corollary IB.5I □ 
By putting Corollary IB.5I and (|A.8|) together we obtain 
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Corollary B.6. With {fj}j=o,i,... defined as in Appendix UQ we have supj ||/j||oo < 00 
and sup.,- \\f"\\oo/j < 00. 

Proof. From (|A.8|) . see also the discussion right after that, we see that /j = —£j/q and 
f" = —ej/q + £jq'/q 2 . Where ej is the j th (normalized) eigenvector. Taking into account 
the normalization, see proof of Corollary IB,5[ the claim is readily proven. □ 

Proof of Proposition \B^j Injecting (IB. lip in the integral term of (lB.16j) leads to 



Vi(0) 



V2p 



le 



V2pi9 



le 



me^e^^dO' 



+ 



(B.41) 



Similarly, we obtain 
y[(6) = y/2pie^ e + 



+ e 



^ e 'me>(e^)de' 



+ o\- 

P 



(B.42) 



and similar expressions for y 2 and y 2 , which actually are just the complex conjugate of y\ 
and y[. We define 



#1 



/*2"7T /*27T 

/ e-^'meie^^dO' , H 2 = / e^ 2 * 6 ' m e > i/ 2 ^) dO' and Q = e 2 ^ 2 ^ 
Jo Jo 



With the higher estimates (|B.41|) and (|B.42|) . we see that (|B.28P becomes 



n-i-^- p [mHx + itiH 2 } + o J -o + 1 - ^ [inHx + inS 2 ] + o 

which implies 



(B.43) 

0, 
(B.44) 



\Q-1\ 2 -^^({n-l)^) = O (J^j . (B.45) 

We now use the expansion of p given by (|B.3ip . In particular, the 0(1/ p 2 ) above becomes 
a 0(l/k 2 ). The second term of the left hand side above is of order 1/k 2 . In fact, we get 
the first order of Hi : 

,-2tx 



c ki0 'me'(e ki -)de' + ( 



where the non local terms in the integral are negligible, since we have 



iK 



^(e ie I k ^(Kr)-e-^I k+1 (Kr)) 



(B.46) 



(B.47) 



2(27rIo(2 J FTr)) 1 /2 

and we can apply (IB.37j) . A similar bound apply for J' * (\fqo~e kt '). So it remains the (real 
!!) first order (remark that J * qo(-) = — Kr sin(-)): 



Hi= JJ^((J*qo) 2 + J*q )(0')d9' + o(^) 



+ 



1 



VkJ ' 



(B.48) 
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But since (using (|B.3ip ) 



n - 1 = 2m(V2f> -k) + 



(B.49) 



where the first term of the right hand side is of order 1/y/k, we have improved the result 
of Lemma EU since using ([R45]) . (|R48jl . (|R49]) and ([BTST]) we obtain 



which implies 



^ = *+°(i) 

Taking (|B.5ip into account, (|B.44j) yields 

2 1 



n- i\ z - -9i((n - + ^(l^ir - 1^1 



(B.50) 



(B.51) 



(B.52) 



The non local terms in H 2 are negligible as for H\ (see above) and a direct calculation 
shows that the local terms are of order 1/k, so from (|B.52p . (|B.48|) and (|B.49P we get 



K 2 r 2 K A r A ■ - 

(V2p - kf ^—(V^P -k) + -—^ = ( V2p - k 



2k 



IQk 2 



K 2 r 2 ^ 2 



4k 



O 



k ? > 



which implies 



K 2 r 2 x / 1 \ 

^ = k + —- k + °{W 2 )- 



, (B.53) 



(B.54) 



We now go further in the expansion to prove that the 0(l/k 3 / 2 ) in (EToip is in fact a 
0(1/ k). Using (|B.17|h we get the the second order expansion of y\ (recall (|B.15P and 

yi (2vr) = n+j^ A (2tt, 6 U f )69 1 +J^ M^uM^u 02, fo))60 1 69 2 +o(J-) . 

'" '' (B.55) 

From (|B.5ip . we deduce 



2- 



A o (e 1 ,e 2 jo)d0 2 



1 



V2pie± 



ie 



V2p 

-V2pi9i 



1 e- V2 P ie2 mg 2 (e V2pi -)d9 2 
Jo 



01 



,V2pid 



2 me 2 (e V2 P ,l -)d9 2 



m e2 (e ki -)66 2 -e- kiyi I e kl ^m e2 {e kv )69 



,ki0 2 „ 



ki- 



and since the non local terms are negligible (see (|B.47P ). we get 



-2- 



A (9 1 ,9 2 ,f )69 2 



iKre 



2k 



Kr 

sin 9i + — sin(26»i' 



Kr 



(B.56) 



(B.57) 
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We deduce the following expansion for the third term of the right hand side of ()B.55[) : 
A (27T,9 1 ,A (9 1 ,e 2 ,f ))de 1 d9 2 = 



Jo 



Kr 
Zk 2 



2jt 



I e- ki01 m dl 



i.: i Kr , . Kr 

e kl ( sin • + — sm(2-) - — 



2tt 



n I e ki ^m ei 



e ( sm- + — sm(2-) - — 



Mi\+0 



k 3 



Using similar arguments as before, we get to 



2tt 



Kr 



e- km m 01 ( e kt -( sin(-) + — sin(2-) ) ) d9 1 = O 



2n 



Kr 



e kMl m 01 ( e kl -{ sin(-) + — sin(2-) ) ) d0i = O 



1 

l 



(B.58) 



(B.59) 



(B.60) 



Moreover, the non local terms of mg 1 (e kl '-) are of order 1/k. In fact, these non local terms 
are finite sums of the form 



2tt 



e nM ^{9)9d9, 



(B.61) 



where \m\ is included in [k — 1, fc+1], and it is easy to see that since the Fourier coefficients 
of ^/qo decay very quickly (see (|B.37|) ). (|B.61|) is of order 1/k. So (|B.58|) becomes 



2tt i-2-k 



A (2ir,9 1 ,A (9 1 ,9 2 ,f ))d9 1 d ° 2 



10 JO 

and we deduce from (|B.55|) 



yi(27r) - in(0) = n - 1 - -(Offi - dH 2 ) 



_kW„ (\ 

8k 2 
K A r 4 7r 2 



Similarly, we obtain 
^(2^-^(0) 
V2pi 



n-i- ~(nHi + inH 2 ) 

k 



8k 2 

K 2 r 2 7r 2 
8k 2 



n + o 



l 
l 



(B.62) 
(B.63) 
(B.64) 



Using these new estimates, (|B.44p becomes 



V2p 



O + O 



K*rV 



o - 1 - -jL- [tfiH! + io# 2 ] -fi + 1 - -k- i inH i + ^H 2 ] 



(B.65) 



which leads to 
x/2 



n - 1| 2 - x± 9f((n - 1)^) + ^(|^| 2 - |# 2 | 2 ) 1 ^ Vvr2 



'4 - 20 - 201 



4 2 



A: 1 



The last term of ()B.66p is of order 1 //c 4 since using (|B.5ip we get 

= 1 + i2ir(V2p -k) + (J^j . 



(B.67) 
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Moreover using (|B.5ip we have 

I f-2lT 



H x = JJ e- kie m e (e ki -) d6 + i(V2p - k) ^ - J J e kit W„(e K ') d6 

+ f 2 J e ki6 m (e ki --) + O (j^J . (B.68) 

As before, the non local terms of ms(e fe '-) are of order 1/k, so the last two integrals in 
(|B.68|) are equal up to a correction of order 1/k, and thus (recall (|B.48|) for the first order 
term), using (|B.51|) . 

irK 2 r 2 „ / 1 



Hi = — + O^J. (B-69) 

We deduce that the first term of the second row of (IB.66j) is of order 1/k 4 , and that, using 
(IB37D . 

^-9((fi - l)Hi) = (V2p - k) + O ( ¥ J , (B.70) 



and 



Since |-£f2| is of order 1/k and that ()B.67p implies 



|0- lj 2 = 4tt 2 (v / 2-p) 2 + 0(^) . (B.72) 



(|B.66P becomes 

'2„2 7^4^,4 



and we deduce 



K 2 r 2 K A r 4: ( 1 \ 

(^-*) 2 -^(^-*) + w = °(f)- (R73 > 

K 2 2 i / i \ 

^ = * + — S (BJ4) 

Now we are able to get a second expansion of the eigenvectors: using (|B.4ip . (|B.57p and 
(|B.74p . we get the following expansion for y\ 

(Kri K 2 r 2 i \ / 1 \ 

1 + — sin(0) + sm(2*)] + O J (B.75) 

and 2/2 is the complex conjugate. So if we define w\ and u>2 the real and imaginary parts, 
we get 

Wl (6) = CO s{k6)-^^(^sm9+^sm{2e^+o(^\ , (B.76) 
w 2 (9) = s in(ke) + C ^^(^ S me + ^sm(2e))+o(^) . (B.77) 



Therefore the proof of Proposition IB.4I is complete. □ 
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